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Transverse momentum distributions and yields for tt*, , p and p in p + p collisions at -^i = 200 
and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy 
Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified 
particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter Tinv, mean 
transverse momentum (pr) and yield per unit rapidity dN/dy at each energy, and compare them 
to other measurements at different ^/s in p + p and p + p collisions. We also present the scaling 
properties such as mr scaling, xt scaling on the pr spectra between different energies. To discuss the 
mechanism of the particle production in p+p collisions, the measured spectra are compared to next- 
to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

Single particle spectra of identified hadrons in high en- 
ergy elementary collisions have attracted physicists for 
many decades due to their fundamental nature and sim- 
plicity. Particle production, in general, can be catego- 
rized into two different regimes depending on the trans- 
verse momentum of the hadrons. One is soft multipar- 
ticle production, dominant at low transverse momentum 
(pT < 2 GeV/c), which corresponds to the ~ 1 fm scale of 
the nucleon radius described by constituent quarks. An- 
other regime is hard-scattering particle production, evi- 
dent at high transverse momentum [p^ > 2 GeV/c) due 
to the hard scattering of point-like current quarks, which 
corresponds to a very short distance scale ^ 0.1 fm [1] 
and contributes less than a few percent of the cross sec- 
tion for y/s < 200 GeV. These two different regimes of 
particle production in p -\- p collisions indicate that "ele- 
mentary" p-\-p collisions are actually rather complicated 
processes. It is interesting to know where the "soft-hard 
transition" happens, and its beam energy and particle 
species dependences, since they have not yet been fully 
understood. 

In soft particle production, cosmic ray physicists ob- 
served in the 1950s that the average transverse momen- 
tum of secondary particles is limited to ~0.5 GeV/c, in- 
dependent of the primary energy Cocconi, Koester 
and Perkins Q then proposed the prescient empirical for- 
mula for the transverse momentum spectrum of meson 
production: 

-^^Ae-^PT, (1) 
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where pT is the transverse momentum in GeV/c and 
(pj,) = 2/6 =0.333 GeV/c. The observation by Orear [| 
that large angle p + p elastic scattering measurements at 
AGS energies (10 to 30 GeV in incident energy) "can be 
fit by a single exponential in transverse momentum, and 
that this exponential is the very same exponential that 
describes the transverse momentum distribution of pions 
produced in nucleon-nucleon collisions" , led to the inter- 
pretation that particle production is "statistical" with 
Eq.[T]as a thermal Boltzmann spectrum, with 1/6=0.167 
GeV/c representing the "temperature" T at which the 
mesons or protons are emitted Q- 

It was natural in a thermal scenario [1, Q to represent 
the invariant cross section as a function of the rapidity 
(y) and the transverse mass (m^ = \/pt^ + 'ti^) with 
a universal temperature parameter T . This description 
explained well the observed successively increasing (pr) 
of TT, K, p, A with increasing rest mass [T^-[l3|, and had 
the added advantage of explaining, by the simple factor 
g~6(mK-m,) _ Y2% , the low value of -10% observed 
for the Kj-K ratio at low pT at ISR energies {^fs ^-^20-60 
GeV) [H. 

In 1964, the constituent quark model with SU(3) sym- 
metry was introduced to explain the hadron flavor spec- 
trum and the static properties of hadrons [3, [l^ . Later 
on, a dynamical model was developed to calculate the 
flavor dependence of identified hadrons in soft multi- 
particle production \\M , together with the inclusive reac- 
tion formalism |17H19| . These theoretical studies on the 
particle production mechanism showed that there was 
much to be learned by simply measuring a single particle 
spectrum, and it brought the study of identified inclusive 
single particle production into the mainstream of p -I- p 
physics. 

One of the controversial issues in understanding soft 
multi-particle production in the 1950's was whether more 
than one meson could be produced in a single nucleon- 
nucleon collision ("multiple production"), or whether the 
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multiple meson production observed in nucleon-nucleus 
{p + A) interactions was the result of several successive 
nucleon-nucleon collisions with each collision producing 
only a single meson ("plural production") 20]. The is- 
sue was decided when multiple meson production was 
first observed in 1954 at the Brookhaven Cosmotron in 
collisions between neutrons with energies up to 2.2 GeV 
and protons in a hydrogen filled cloud chamber [d, f2l| . 

Then the observation of multi-particle production oc- 
curring not only in nucleon-nucleus ij) -\- A) but also in 
nucleon-nucleon {p -t- p) collisions motivated Fermi and 
Landau to develop the statistical [22| and hydrodynam- 
ical j23| approach to multi-particle production. Belenkij 
and Landau observed that although the statistical model 
of Fermi is sufficient to describe the particle numbers 
in terms of only a temperature and a chemical poten- 
tial, this model has to be extended to hydrodynamics, 
when particle spectra are considered. They also noted 
that the domain of the applicability of ideal relativistic 
hydrodynamics coincides with the domain of the appli- 
cability of thermodynamical models in high energy p -\- p 
collisions [23]. 

On the other hand, understanding of the particle pro- 
duction by hard scattering partons has been advanced 
by the appearance of a rich body of data in p -\- p colli- 
sions at the CERN-ISR [Hill El in the 1970s, followed 
by measurements at the Relativistic Heavy Ion Collider 
(RHIC) dXyfs = 200 [26 33] and 62.4 GeV H over the 
last decade. The hard scattering \n p + p collisions has 
been discovered by the observation of an unexpectedly 
large yield of particles with large transverse momentum 
and the phenomena of dijets at ISR [ssj. These obser- 
vations indicate that hard scattering process occurs be- 
tween the quarks and gluons constituents (or partons) 
inside the nucleons. This scattering process can be de- 
scribed by the perturbative quantum chromodynamics 
(pQCD) because the strong coupling constant as of QCD 
becomes small (asymptotically free) for large momentum 
transfer (Q^) parton-parton scatterings. After the initial 
high parton-parton scatterings, these partons frag- 
ment into the high pt hadrons or jets. In fact, at RHIC 
energies, single particle spectra of high pT hadrons are 
weU described by pQCD [s^, [H H]. Furthermore, xt 
(= 2pTl\fs), which is also inspired by pQCD, is known 
to be a good scaling variable of the particle production 
at high PT at both ISR ^ and RHIC [sS] energies, so 
that XT scaling can be used to distinguish between the 
soft and hard particle productions. 

Another important point of measurements in p~\-p colli- 
sions is as a baseline for the heavy ion (A+A) data. The 
nuclear modification factor i?AA, for example, uses pt 
spectra in p -f p collisions as a denominator, and those in 
A-l-A collisions (with the appropriate scaling of number 
of binary nucleon-nucleon collisions) as a numerator. In 
addition, pT spectra in p -I- p provide a reference of bulk 
properties of A-l-A collisions, such as the inverse slope 
parameter Tinv, mean transverse momentum (pr), and 
yield per unit rapidity dN / dy. These data in p + p col- 



lisions can be treated as baseline values for the smallest 
A+A collisions. 

In this paper, we present measurements of identified 
charged hadron pT spectra for tt^, , p and p at midra- 
pidity in p + p collisions at -^s = 200 and 62.4 GeV from 
the PHENIX experiment. First, we compare the results 
of particle spectra at 200 GeV with those at 62.4 GeV 
as a function of pt, "ITj and tot — "m (where m is the 
rest mass). Second, the extracted values from pt spec- 
tra, i.e. Tinv, (pt), and dN/dy, are compared between 
the two beam energies. For the systematic study of par- 
ticle production as a function of the data are further 
compared to measurements in p + p and p -\'P collisions 
at the CERN-ISR and FNAL-Tevatron colliders. 

From these measurements, we discuss the following key 
issues: 

hard scattering particle production - The data are 
compared with the results of perturbative quantum 
chromodynamics (pQCD) calculations. 

transition from soft to hard physics - Since the pt 
regions presented in this paper can cover the region 
where the soft-hard transition occurs, the scaling 
properties in ttit and xt with their beam energy 
and particle species dependences are shown. 

comparisons with heavy-ion data as a baseline 
measurement - Some of the data in p + p are 
compared with the existing data in Au+Au 37]. 

The paper is organized as follows: Section HIl describes 
the PHENIX detector as it was used in this measure- 
ment. Section lHll discusses the analysis details, including 
data sets, event selection, track selection, particle identi- 
fication, corrections applied to the data, and systematic 
uncertainties. Section IIVI gives the experimental results 
for Pt spectra for identified charged particles, particle 
ratios, rriT scaling, the excitation function of observables 
(such as Tinv, (pt), and dN/dy), and -Raa- Section IVl 
compares the results with next-to-leading-order (NLO) 
and next-to-leading-log (NLL) pQCD calculations, and 
discuss soft and hard particle production, and the tran- 
sition between them. Section [VII gives the summary and 
conclusions. 



II. EXPERIMENTAL SETUP 

The PHENIX experiment is designed to perform a 
broad study of A-fA, d + A, a.nAp + p collisions to inves- 
tigate nuclear matter under extreme conditions, as well 
as to measure the spin structure of the nucleon. It is 
composed of two central arms (called east and west arm, 
respectively), two forward muon arms, and global de- 
tectors, as shown in Fig. [TJ The central arms are de- 
signed to detect electrons, photons and charged hadrons 
in the pseudorapidity range |?7| < 0.35. The global detec- 
tors measure the start time, collision vertex, and charged 
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West Beam View East 




FIG. 1: (color online) The PHENIX detector configuration 
for RHIC Run-6 data taking period. The beam-beam counter 
(BBC) is labeled as BB on the figure. 



hadron multiplicity of the interactions in the forward 
pseudorapidity region. The following sections describe 
those parts of the detector that are used in the present 
analysis. A detailed description of the complete set of 
detectors can be found elsewhere (38l - l42| . 

The beam-beam counters (BBC) [4l| determine the 
start time information for time-of-flight measurements 
and the collision vertex point, as well as providing the 
main collision trigger. The two BBCs are located at 1.44 
m from the nominal interaction point along the beam 
line on each side. Each BBC comprises 64 Cerenkov 
telescopes, arranged radially around the beam line. The 
BBCs measure the number of charged particles in the 
pseudorapidity region 3.0 < \ri\ < 3.9. 

Charged particle tracks are reconstructed using the 
central arm spectrometers 42]. The east arm spectrom- 
eter of the PHENIX detector contains the following sub- 
systems used in this analysis: drift chamber (DC), pad 
chamber (PC) and time of flight (TOF). The magnetic 
field for the central arm spectrometers is supplied by the 
central magnet [l^ that provides an axial field parallel 
to the beam around the collision vertex. 

The drift chambers are the closest tracking detectors 
to the beam line, located at a radial distance of 2.2 m (ge- 



ometric center, same for other detectors). They measure 
charged particle trajectories in the azimuthal direction to 
determine the transverse momentum of each particle. By 
combining the polar angle information from the first layer 
of PC, as described below, with the transverse momen- 
tum, the total momentum p is determined. The momen- 
tum resolution in p+p collisions is 6p/p ~ 0.7%®1.0%xp 
(GeV/c), where the first term is due to the multiple scat- 
tering before DC and the second term is the angular reso- 
lution of the DC. The absolute momentum scale is known 
as ± 0.7% rms from the reconstructed proton mass using 
TOF. 

The pad chambers are multi-wire proportional cham- 
bers that form three separate layers of the central track- 
ing system. The first layer (PCI) is located at the radial 
outer edge of each drift chamber at a distance of 2.49 m, 
while the third layer is at 4.98 m from the interaction 
point. The second layer is located at a radial distance of 
4.19 m in the west arm only. PCI and DC, along with 
the vertex position measured by BBC, are used in the 
global track reconstruction to determine the polar angle 
of each charged track. 

The time-of-flight detector serves as the primary parti- 
cle identification device for charged hadrons by measur- 
ing the stop time. The start time is given by BBC. TOF 
is located at a radial distance of 5.06 m from the inter- 
action point in the east central arm. This contains 960 
scintillator slats oriented along the azimuthal direction. 
It is designed to cover I77I < 0.35 and Acj) = 45° in az- 
imuthal angle. The intrinsic timing resolution is cr ~ 115 
ps, which in combination with the BBC timing resolu- 
tion of 60 ps, allows for a 2.6cr n/K separation at px — 
2.5 GeV/c, and K/p separation out to pr = 4.5 GeV/c, 
using an asymmetric particle- identification (PID) cut, as 
described below. 



III. DATA ANALYSIS 

The two RHIC data sets analyzed are 2005 data for 
p + p coUisions at ^/s = 200 GeV and 2006 data for 
p + p collisions at ^/s — 62.4 GeV. Each data set was 
analyzed separately by taking into account the different 
run conditions and accelerator performance. In this sec- 
tion, we explain the event selection, track reconstruction, 
particle identification, and corrections to obtain the px 
spectra. The event normalization and systematic uncer- 
tainties are also presented. 



A. Event selection 

We use the PHENIX minimum bias trigger events, 
which arc determined by a coincidence between north 
and south BBC signals, requiring at least one hit on both 
sides of BBCs. Due to the limited acceptance, approxi- 
mately only half of the p + p inelastic cross section can 
be measured by BBC. The PHENIX minimum bias data. 
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triggered by BBC in p+p collisions within a vertex cut of 
±30 cm, include ctbbc = 23.0 ±2.2 nib at = 200 GeV 
and ctbbc — 13-7 ±1.5 mb at ^/s — 62.4 GeV (see Sec- 
tion |IITE|. We analyze 9.2x10^ minimum bias events for 
the 2005 p+p data at -^/s — 200 GeV, which is more than 
30 times larger than the 2003 data set [111 and 2.14x10* 
times larger than the minimum bias events for the 2006 
data at = 62.4 GeV. 



B. Track reconstruction and particle identification 

As in previous publications [13, ' charged particle 
tracks are reconstructed by DC based on a combinatorial 
Hough transform, which gives the angle of the track in 
the main bend plane. PCI is used to measure the po- 
sition of the hit in the longitudinal direction along the 
beam axis. When combined with the location of the col- 
lision vertex along the beam axis, the PCI hit gives the 
polar angle of the track. Only tracks with valid infor- 
mation from both DC and PCI are used in the analysis. 
To associate a track with a hit on TOF, the track is pro- 
jected to its expected hit location on TOF. We require 
tracks to have a hit on TOF within ±2(t of the expected 
hit location in both the azimuthal and beam directions. 
The track reconstruction efficiency is approximately 98% 
in p + p collisions. Finally, a cut on the energy loss in 
the TOF scintillator is applied to each track. This f3- 
dependent energy loss cut is based on a parameterization 
of the Bethe-Bloch formula. The flight path length is cal- 
culated from a fit to the reconstructed track trajectory 
in the magnetic field. The background due to random 
association of DC/PCl tracks with TOF hits is reduced 
to a negligible level when the mass cut used for particle 
identification is applied. 

Charged particles are identified using the combination 
of three measurements: time-of-flight from the BBC and 
TOF, momentum from the DC, and flight-path length 
from the collision vertex point to the TOF hit position. 
The square of mass is derived from 



P 



/ itof 



1 



(2) 



where p is the momentum, itof is the time of flight, L is 
the flight path length, and c is the speed of light. The 
charged particle identification is performed using cuts in 
and momentum space. In Fig. ^ a plot of momen- 
tum multiplied by charge versus is shown together 
with applied PID cuts as solid curves. We use 2cr stan- 
dard deviation PID cuts in and momentum space for 
each particle species. The PID cut is based on a param- 
eterization of the measured width as a function of 
momentum 



^(4mV) 



£2 



[4p2 



Am 



where ctq is the angular resolution, ams is the multiple 
scattering term, at is the overall time-of-flight resolu- 
tion, TO is the centroid of distribution for each par- 
ticle species, and Ki is the magnetic field integral con- 
stant term of 101 mrad-GeV. The parameters for PID are 
(Tq = 0.99 mrad, ams — 1-02 mrad-GeV, and at — 130 ps. 
For pion identification above 2 GeV/c, we apply an asym- 
metric PID cut to reduce kaon contamination of pions. 
As shown by the lines in Fig. [2l the overlap regions which 
are within the 2a cuts for both pions and kaons are ex- 
cluded. The lower momentum cut-offs are 0.3 GeV/c for 
pions, 0.4 GeV/c for kaons, and 0.5 GeV/c for protons 
and antiprotons. The lower momentum cut-off value for 
p and p is larger than for pions and kaons due to the 
larger energy loss effect. 

For kaons, the upper momentum cut-off is 2 GeV/c 
since the tt + p contamination level for kaons is 8% 
at that momentum. The upper momentum cut-off for 
pions is pt — 3 GeV /c where the K + p contamination 
reaches ~ 3%. Electron (positron) and decay muon back- 
ground at very low pT (< 0.3 GeV/c) are well separated 
from the pion mass-squared peak. For protons the up- 
per momentum cut-off is set at 4.5 GeV/c. For protons 
and antiprotons an additional cut, > 0.6{GeV /c^)^ , 
is introduced to reduce the contamination. The contam- 
ination background on each particle species is subtracted 
statistically after applying these PID cuts. 
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FIG. 2: (color online) Momentum multiplied by charge versus 
mass squared distribution in p+p collisions at ^/s = 62.4 GeV. 
The lines indicate the PID cut boundaries {2a) for pions, 
kaons, and protons (antiprotons) from left to right, respec- 
tively. 



C. Efficiency corrections 

We use the PHENIX-integrated-simulation applica- 
tion, which is a GEANT [44] based Monte Carlo (MC) 
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simulation program of the PHENIX detector, to correct 
for geometrical acceptance, reconstruction efficiency, in- 
flight decay for tt and multiple scattering effect, and 
nuclear interactions with materials in the detector (in- 
cluding p absorption). Single particle tracks are passed 
from GEANT to the PHENIX event reconstruction soft- 
ware [43i]- In this simulation, the BBC, DC, and TOE 
detector responses are tuned to match the real data. Eor 
example, dead areas of DC and TOE are included, and 
momentum and time-of-ffight resolutions are tuned. The 
track association to TOE in both azimuth and along the 
beam axis as a function of momentum and the PID cut 
boundaries are parameterized to match the real data. A 
flducial cut is applied to choose identical active areas on 
TOE in both the simulation and data. 

We generate 1x10^ single particle events for each par- 
ticle species {t:^ , , p and p) with flat pr distributions 
for high pt (2-4 GeV/c for pious and kaons, 2-8 GeV/c 
for p and p) with enhancement at low pT (< 2 GeV/c). 
Weighting functions to the pr distributions are also used 
to check the effect of steepness, which is less than ^1% 
level on the flnal yields in the measured pt range. The 
rapidity range is set to be wider than the PHENIX ac- 
ceptance, i.e. flat in -0.6 < y < 0.6 (Ay = 1.2) to deal 
with particles coming from outside (the denominator of 
Eq. [4]is weighted with a factor l/Ay = 1/1.2 in order to 
normalize the yield for unit rapidity). The efficiencies are 
determined in each pT bin by dividing the reconstructed 
output by the generated input as expressed as follows: 



3. measured p jp) distributions scaled with ISR A (A) 
distributions [25|. 

Using each input above, proton and antiproton spectra 
from weak decays are calculated by using Monte Carlo 
simulation to take into account decay kinematics, the 
PHENIX track reconstruction efficiency and experimen- 
tal acceptance. Then systematic uncertainties are evalu- 
ated from different A and A spectra inputs. The resulting 
uncertainties on the flnal proton and antiproton spectra 
are of the order of 20-30% at pr = 0.6 GeV/c and 2- 
5% at = 4 GeV/c. The fractional contribution of the 
feed-down protons (antiprotons) to the total measured 
proton (antiproton) spectra, (5focd(PT), is approximately 
10-20 (5-15) % at PT = 4 GeV/c for 200 GeV p+p (62.4 
GeV p+p) and it shows an increase at lower pT as shown 
in Eig. [3l The correction factor for the feed-down correc- 
tion can be expressed as: Ctccd(PT) = 1 — <5fcod(PT), which 
is multiplied to the raw pT spectra (see Section IIII Gp . 

The feed-down correction for protons is different from 
that for antiprotons in 62.4 GeV, because of the difference 
in h./p and h/p ratio at this beam energy. At 62.4 GeV, 
A/p ratio is 0.2 while A/p ratio is 0.4 [2^, so that 
the feed-down contribution for antiprotons is bigger than 
that for protons. At 200 GeV, these two ratios are almost 
same [3l|, therefore the feed-down corrections for p and 
p become identical. 

E. Cross section normalization 



# of reconstructed MC tracks 
^^^^-^ " # of generated MC tracks ' ^ ' 

The resulting correction factors Coff(pT) (= l/e(pT)) 
is multiplied to the raw px spectra for each px bin and 
for each individual particle species (see Section UlI G[) . 



D. Feed-down corrections 

The proton and antiproton pT spectra are corrected 
for feed-down from weak decays of hyperons. The de- 
tailed procedure of the feed-down correction can be found 
in [2^. We include the following decay modes: A p7r~, 
S+ — )■ p7r°, and A production from E", S~. The 
feed-down contributions for antiproton yields are also es- 
timated using the above decay modes for antiparticles. 

In order to estimate the fractions of protons and an- 
tiprotons from weak decays of hyperons in the measured 
proton and antiproton pT spectra, we use three input A 
and A pr spectra: 

1. measured A and A pT spectra in PHENIX in p + p 
collisions at = 200 and 62.4 GeV, 

2. measured p (p) distributions scaled with measured 
A (A) distributions ^31i] . and 



The BBC counter serves a dual function as both the 
minimum bias trigger and the calibrated luminosity mon- 
itor. The luminosity C is deflned as the interaction rate 
for a given cross section: dN/dt = Ca and the total num- 
ber of events for a given cross section is 

N = ax J Cdt, (5) 

where J Cdt is the integrated luminosity. To connect 
the number of minimum bias triggered events and the 
integrated luminosity, ctbbc is introduced where 1 / ctbbc 
corresponds to the integrated luminosity per minimum 
bias triggered event (Eq. IH]). 

A^BBC = CTBBC X J Cdt, (6) 

where iVBBC is the number of minimum bias events and 
J Cdt is the corresponding integrated luminosity. (Tbbc 
is measured by a Van der Meer scan method (Vernier 
scan) in PHENIX [13, SI. 

Vernier scans were performed for ^/s — 200 and 62.4 
GeV data sets. The obtained aBBC are 23.0 ± 2.2 mb 
and 13.7 ± 1.5 mb for ^ = 200 and 62.4 GeV respec- 
tively. The quoted uncertainty is a systematic uncer- 
tainty. These numbers were reported in our measure- 
ments of 7r° production [s^, [131 ■ 
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FIG. 3: (color online) Fraction of feed-down protons and an- 
tiprotons as a function of pr with systematic uncertainties. 
Top: 200 GeV p-l-p (positive and negative functions are com- 
mon). Bottom: 62.4 GeVp + p. 



trigger bias. We use this pr dependent trigger bias cor- 
rection for charged hadrons, by using fitted coefficients 
of a second order polynomial, as shown in Fig. |4l 

With those values, the invariant yield per BBC trigger 
count {Y/Nqbg) is related to the invariant cross section 
(ct) using 
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FIG. 4: (color online) Fraction of the inclusive vr'' yield which 
satisfied the BBC trigger condition in 62.4 GeV p + p- Data 
points are from Fig. 1 of [S^ ]. 



F. Systematic uncertainties 



Since the minimum bias trigger registers only half of 
the p + p inelastic cross section, it is expected that there 
is a trigger bias against particles in the central spectrom- 
eters. This was checked with tt'^ 's in the electromagnetic 
calorimeter with high px photon triggered events, and 
with charged tracks in accelerator's beam crossing (clock) 
triggered events. The trigger bias ebias determined from 
the ratio (/^o) of the number of tt^ in the high px pho- 
ton triggered sample with and without the BBC trigger 
requirement [S^l . We assume ebias is process dependent 
and so that it is measured as ebias = /tt"- This ratio, 
/^o, is 0.79 ±0.02 independent of the transverse momen- 
tum for = 200 GeV. At 62.4 GeV, the trigger bias 
was found to be transverse momentum dependent [s^ . 
Fig. m shows that the trigger bias Jt^o, and Jt^o is «40% 
up to pt ~3 GeV/c, and monotonically decreases to 25% 
at ~7 GeV/c. As described in the previous PHENIX 
publication [3J|, this decrease can be understood by the 
fact that most of the energy is used for the production of 
high-energy jets which contain the measured high pT tt'^ 
and charged hadrons, and there is not enough energy left 
to produce particles for y/s = 62.4 GeV p+p collisions at 
the forward rapidity (3.0 < \ri\ < 3.9) where the BBC is 
located. This drop can be seen only for 62.4 GeV data. 
Also we assume no particle species dependence for this 



In order to estimate the systematic uncertainties, the 
Pt spectra with slightly different analysis cuts from what 
we use for the final results are prepared, and these spectra 
are compared to those with the standard analysis cuts. 
We checked the following analysis cuts: (1) fiducial, (2) 
track association windows, and (3) PID. 

For each spectrum with modified cuts, the same 
changes in the cuts are made in the Monte Carlo sim- 
ulation. The fully corrected spectra with different cut 
conditions are divided by the spectra with the baseline 
cut condition, resulting in uncertainties associated with 
each cut condition as a function of pr- The obtained 
uncertainties are added in quadrature. Tables U and [Tl] 
show the systematic uncertainties on px spectra for each 
data set. There are three categories of systematic un- 
certainty: Type A is a point-to-point error uncorrelated 
between px bins, type B is px correlated, where all points 
move in the same direction but not by the same factor, 
while in type C all points move by the same factor inde- 
pendent of Pt [46j . In this study, the systematic uncer- 
tainties on feed-down correction and PID contamination 
correction are Type B, other systematic uncertainties on 
applied analysis cuts are Type C. There are two types 
of the PID-related uncertainties. One is the systematic 
uncertainty of the yield extraction, which is evaluated by 
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TABLE I: Systematic uncertainties on the pr spectra for 
^/s = 200 GeV p + p given in percent. The number in paren- 
thesis includes the pr dependence of the uncertainties for PID 
cut, feed-down correction, PID contamination correction. 



TABLE IL Systematic uncertainties on the pr spectra for 
^/s = 62.4 GeV p + p given in percent. The number in paren- 
thesis includes the pr dependence of the uncertainties for feed- 
down correction, PID contamination correction. 



Source 


ir+ 


TV 


K+ 


K- 


P 


P 


Source 


7r+ 


TT 


K+ 




P 




P 


Fiducial Cut 


5 


5 


4 


5 


4 


5 


Fiducial Cut 


6 


5 


6 


5 


7 




5 


Track Matching 


4 


4 


5 


4 


4 


4 


Track Matching 


2 


2 


3 


3 


3 




3 


PID Cut 


3 


3 


2 


2 


2-8 


2-10 


PID Cut 


2 


2 


3 


3 


4 




4 


Efficiency Correction 


2 


2 


2 


2 


2 


2 


Efficiency Correction 


2 


2 


2 


2 


2 




2 


Feed-down Correction 










4-25 


4-25 


Feed-down Correction 










1- 


-16 


3-50 


PID Contamination 










0-2 


0-2 


PID Contamination 






0-5 


0-5 








Total 


7 


7 


7 


7 


6 (8-25) 


7 (9-25) 


Total 


7 


6 


7 


7 


9 


(9-18) 


7 (8-50) 



changing the PID boundary in vs monientum plane. 
The other is the systematic uncertainty of the particle 
contamination, which is evaluated by using the contam- 
ination fraction. The fraction is estimated by fitting m? 
distributions on each pT slice under the several condi- 
tions as follows: (1) fixed parameters for p and p mass 
centroid and width, (2) p and p mass centroid free with 
fixed mass width, (3) p and p mass width free with fixed 
mass centroid. 

The systematic uncertainty on the BBC cross section 
is 9.7% and 11% for = 200 and 62.4 GeV respectively. 
The systematic uncertainty on the trigger bias is 3% and 
1-5 % for y/s =^ 200 and 62.4 GeV respectively (see Sec- 
tion |IITE]). These uncertainties on normalization (Type 
C) are not included in Tables HI and HIl All the figures and 
tables including the tables in Appendix [Xl do not include 
the normalization uncertainties, unless explicitly noted. 



G. Invariant cross section 



The differential invariant cross section is determined 



E 



, d^cr _ ^BBC 

dp3 ~ 2^pTiVBBcCgBC(p^) 



Coff(pT)Cfocd(PT) 



dprdy' 
(8) 

where a is the cross section, pT is transverse momentum, 
and y is rapidity, A^bbc is the number of minimum bias 
events, ubbc is the minimum bias cross section measured 
by BBC, Ccs{pt) is the acceptance correction factor in- 
cluding detector efficiency, C^i^ipr) is the trigger bias, 
CiccdiPr) is the feed-down correction factor only for pro- 
tons and antiprotons, and N is the number of measured 
tracks. 



IV. RESULTS 

In this section, we show the transverse momentum dis- 
tributions and yields for tt^ , , p and p in p -\- p col- 
lisions at = 200 and 62.4 GeV at midrapidity mea- 



sured by the PHENIX experiment. We also present the 
transverse mass (my) spectra, the inverse slope parame- 
ter Tinv, mean transverse momentum (pt), yield per unit 
rapidity dN/dy, and particle ratios at each energy, and 
compare them to other measurements at different -^/s in 
p + p and p + p collisions. The measured Tinv, (pt), and 
dN/dy in p + p 200 GeV are also compared with those 
in published results in Au-|-Au at 200 GeV. The nuclear 
modification factor Raa for 200 GeV Au-|-Au using the 
present study in p + p 200 GeV are also presented. 



A. pt spectra 

Figure [5] shows transverse momentum spectra for tt^ , 
K^, p, and p in 200 and 62.4 GeV p + p collisions. Feed- 
down correction for weak decays is applied for p and p, 
and the same correction factors are consistently used for 
all figures through Section |IV] unless otherwise specified. 
The systematic uncertainty on the BBC cross section is 
9.7% and 11% for = 200 and 62.4 GeV respectively. 
Each of Pt spectra is fitted with an exponential func- 
tional form: 



1 



2ttpt dydpT 



(9) 



where A is a normalization factor and T is an inverse 
slope parameter for px- The fitting parameters and 
xVNDF by using Eq. [S] for 7r±, , p, and p in 200 
and 62.4 GeV p + p collisions, are tabulated in Table HlIl 
The fitting range is fixed as pr = 0.5-1.5 GeV/c for tt^, 
0.6-2.0 GeV/c for K^, and 0.8-2.5 GeV/c for p, p at 
both collision energies. 

Figure [5] shows that pious, protons, and antiprotons 
exhibit an exponential spectral shape at low px and a 
power-law shape at high px, while kaons are exponential 
in the measured px range. The transition from exponen- 
tial to power-law can be better seen at ~ 2 GeV/c 
for pions, and at ~ 3 GeV/c for protons and antipro- 
tons at both energies. The fractions of soft and hard 
components gradually change in the transition region. 
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FIG. 5: (color online) (top, middle) Transverse momentum distributions for tt^, K"^, p, and p in p + p collisions at y/s = (left) 
200 and (right) 62.4 GeV at midrapidity. Only statistical uncertainties are shown, (middle plots) Each spectrum is fitted with 
an exponential function, (lower panels of middle plots) Ratio of the exponential fit to data for each particle species, (bottom) 
Ratios of pt spectra for n"^ , tt" [s^, [1^, K"^, p, and p in 200 GeV p + p collisions to those in 62.4 GeV p + p collisions. 
Statistical and systematic uncertainties are combined in quadrature. The trigger cross section uncertainty is not included. The 
lines represent the NLO pQCD calculations (DSS fragmentation function) for pious with different factorization, fragmentation, 
and renormalization scales (which are equal) [i?! ]. 



Ratios of the px spectra at 200 GeV to those at 62.4 
GeV are shown in the bottom plot of Fig. [S) The left 
panel shows the ratios for positively charged particles, 
and the right panel is those for the negatively charged 



particles. The data for neutral pions [13, [13] are also 
shown on both panels. The ratios show a clear increase 
as a function oipx for all the ratios. Since hard scattering 
is expected to be the dominant particle production pro- 
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FIG. 6: (color online) Transverse mass distributions for tt^, n'^, , p, and p in p + p collisions at = (left) 200 and (right) 
62.4 GeV at midrapidity for (upper) positive and (lower) negative hadrons. Only statistical uncertainties are shown. The 
references for STAR data are 7r='=, p, p [33[] and iC^, A, A [sj]. The references for ISR data are tt*, , p, p and A, A [25l ]. 



cess at high p^, this strong dependence indicates t^wo 
features: (1) spectral shape is harder for 200 GeV com- 
pared to those for 62.4 GeV, and (2) a universal shape 
for all particle species up to pr = 2-3 GeV/c. In the 
same figure, the results from NLO pQCD calculations 
-with DSS fragmentation function [47] for pions "with dif- 
ferent factorization, fragmentation, and renormalization 
scales ("which are equal) are also sho"wn. The agreement 
is relatively poor, due to the disagreement bet"ween the 
NLO pQCD calculation "with DSS fragmentation func- 
tion and measurement for pions at ^/s — 62.4 GeV. As 
"we -will discuss in detail in Section IV C[ it is found that 
NLL pQCD gives a better description of the data for p+p 
at 62.4 GeV. 



Please note that each line in pQCD is calculated for 
each ^ {= pt/2, pT, 2pt) value. The hard scale resides 
in the hard scattering, "which is expected to be the same 
regardless of hadron species. The theoretical uncertainty 
in the ratio of NLO (200 GeV / 62.4 GeV) significantly 
cancels. The same comparison of ratio for NLL cannot 
be made due to the unreliability of resummation in NLL 
pQCD at 200 GeV in the low pr region [471] . 



B. rriT spectra 

In p -I- p{p) collisions at high energies, the transverse 
mass (ttjt) spectra of identified hadrons sho"w a universal 
scaling behavior, and this fact is kno"wn as tot scaling. 
In order to check the mx scaling and to gain a further 
insight on the particle production mechanism especially 
at high pt at RHIC energies, transverse mass spectra in 
200 and 62.4 GeV p+p collisions are sho-wn in Fig. [HI The 
data for tt^, K^, p, and p in 200 and 62.4 GeV are from 
this study. The 7r° spectra are taken from the PHENIX 
measurements [30ll34]. From the STAR experiment, tt^, 
p, and p spectra in 200 GeV p + p are taken from [ssl ]. 
and K^,A, and A spectra in 200 GeV p + p are taken 
from 31|. The tt^, K^, p, p spectra in 63 GeV p+p 
are from [l^, and A, A spectra in 63 GeV p + p are 
from the ISR experiment [25|. For both energies one 
can see "almost" the same spectral shape for all particle 
species. Ho"wever, this scaling is not perfect, because it is 
violated by the yields that are different for each particle 
species. In order to see the approximate scaling behavior 
and its violation more clearly, one needs to introduce 
certain factors, that normalize these yields to make them 
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FIG. 7: (color online) Scaled transverse mass distributions for tt , tt", K , p, and p in p + p collisions at ^/s — (upper) 200 
and (lower) 62.4 GeV at midrapidity for (upper left) positive, (upper right) negative, and (lower) ± hadrons. Only statistical 
uncertainties are shown, (upper) The STAR spectra for K^, A, A are from 31]. (lower) The ISR spectra for A, A are from 25]. 
Arbitrary scaling factors are applied to match the yield of other particles to that of charged pions in the range of mr ~ 1.0-1.5 
GeV/c^. The lower panels of each plot show the ratio to tt'' Tsallis fit. 



coincide at a certain value of mx- 

Figure [7] shows the ttit spectra with such scaling fac- 
tors implemented. These normalization scaling factors 
are determined to match the yield of each particle species 
to that of charged pions in the range of tot = 1.0-1.5 
GeV/c^, i.e. factors which scale the normalization (verti- 
cal axis) as in the legends of Fig. O These scaling factors 
are tabulated in Table IIVI The bottom panels on the 
plots in Fig. [7] are the ratio of data to the fitting result 
using Tsallis function for 7r° data in 200 GeV [30| and 
62.4 GeV Above my > 1.5 GeV/c^, these figures 

indicate a clear separation between meson and baryon 
spectra. The meson spectra are apparently harder than 
the baryon spectra in this representation. This effect can 
be seen more clearly on the ^/s = 200 GeV dataset, than 
on data measured at 62.4 GeV. Such a baryon-meson 
splitting in tot spectra has been reported by the STAR 
experiment mp+p collisions at y/s = 200 GeV |3l|. That 



only a (quark, antiquark) pair, while baryon production 
requires a (diquark, antidiquark) pair. 



paper argued for a given jet energy, mesons might be pro- 
duced with higher transverse momentum than baryons, 
because meson production in jet fragmentation requires 



However, such a rescaling along the vertical axis is arbi- 
trary. Since it is known that the shapes of the power-law 
tail for different particle species at high pT in 200 GeV 
data are very similar p3 | , one can show by using different 
scaling factors to^ spectra match at higher rriT- In this 
case, rriT spectra for baryons overshoot those for mesons 
at low mx- The interpretation of rriT scaling is not settled 
yet, but from the experimental results presented here, it 
is clear that the spectral shape for mesons and baryons 
are very similar for both low and high rriT regions, but 
different for the intermediate to^ (1.5-2.5 GeV/c^) re- 
gion. In Section HV CI we will discuss the spectral shape 
at low Tot in detail, by taking into account the hadron 
mass effect. 
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TABLE III: Fitting results for A, T of Eq.[9]for pr spectra for 
TT*, K^, p, and p in 200 and 62.4 GeV p + p collisions. The 
fitting range is fixed as pt = 0.5-1.5 GeV/c for tt^, 0.6-2.0 



GeV/c for K"^ 
energies. 



and 0.8-2.5 GeV/c for p, p at both collision 





hadron 


A 


r 


xVndf 




(GeV) 






(GeV/c) 






200 


7r+ 


80.1 ± 7.2 


0.220 ± 0.004 


11.5/8 






n~ 


80.7 ± 7.5 


0.220 ± 0.004 


13.5/8 


S 




K+ 


6.45 ± 0.50 


0.296 ± 0.005 


29.4/12 


d 






6.62 ± 0.51 


0.293 ± 0.004 


18.8/12 


E 




P 


3.24 ± 0.38 


0.318 ± 0.006 


3.3/15 






P 


2.83 ± 0.35 


0.318 ± 0.006 


2.8/15 




62.4 


7r+ 


78.0 ± 7.0 


0.203 ± 0.003 


9.0/8 






7r~ 


81.0 ± 6.2 


0.200 ± 0.003 


11.1/8 


O 




K+ 


6.17 ± 0.52 


0.264 ± 0.004 


15.6/12 


1 




K- 


6.01 ± 0.49 


0.254 ± 0.004 


10.0/12 






P 


4.61 ± 0.48 


0.275 ± 0.005 


2.8/15 


LU 




P 


2.95 ± 0.36 


0.267 ± 0.005 


2.9/15 





TABLE IV: Normalization scaling factors for rriT spectra for 
Fig. [T] The scaling factors for the STAR experiment are de- 
termined from [siT |3^ and those for the ISR results are de- 
termined from [2511. 



Exp. 

(GeV) 



tt" K+ K- Kt 



P 



A A 



200 PHENIX 1.0 1.0 0.9 2.4 2.4 - 1.15 1.4 - - 

200 STAR 1.0 1.0 - - - 2.4 0.75 0.75 0.8 0.9 

62.4 PHENIX 1.0 1.0 0.9 2.32 2.88 - 0.9 1.5 - - 

63 ISR -------- 0.4 0.5 




rriT — m spectra 



FIG. 8: (color online) mT — m spectra for tt^, K^, p, and p in 
p + p collisions at ^/s = (upper) 200 and (lower) 62.4 GeV at 
midrapidity. Only statistical uncertainties are shown. Each 
spectrum is fitted with an exponential form of Eq. (|10|l in the 
range of niT — m = 0.3-1.0 GeV/c^. Solid lines represent the 
functions in the fitted range, dashed lines show the extrap- 
olation of these functions beyond this range, (lower panels) 
Ratio of the exponential fit to data for each particle species. 



Figure |8] shows the mx — m spectra for tt^ , , p, 
and p in 200 and 62.4 GeV p+p collisions, respectively. 
When analyzing these mx — m spectra of various identi- 
fied hadrons, one discusses the spectral shape mainly in 
the low mx — m region. Each of these spectra is fitted 
with an exponential functional form: 

1 cPa . / mx — m\ , , 

TT^^^^^P ^ (10) 

Zirmx ayamx ^ Jinv ^ 

where ^ is a normalization factor and Tinv is called the 
inverse slope parameter. The fitting parameters and 
xVNDF by using Eq. [10] for 7r±, K^, p, and p in 200 
and 62.4 GeV p + p collisions, are tabulated in Table [TTl[ 
The fitting range is fixed as mx — m — 0.3-1.0 GeV/c^ for 
all particle species at both collision energies. We obtain 
smaller x^/NDF for protons and antiprotons than those 
for pious and kaons, because of the larger systematic un- 
certainties for p and p at low px due to the uncertainties 
of weak decay feed-down corrections. As seen in Fig. [5] 



the spectra are exponential in the low mx — m range. At 
higher transverse mass, the spectra become less steep, 
corresponding to an emerging power law behavior. The 
transition from exponential to power law can be seen at 
mx — m = 1-2 GeV/c^ for all particle species. 

The dependence of Tinv on hadron mass is shown in 
Fig. [5| These slope parameters are almost independent of 
the energy of ^/s = 62.4 and 200 GeV. The inverse slope 
parameters of kaons is similar to that of protons while 
the slope parameter of pious has slightly smaller values. 
It may be possible that lower Tjnv values for pious are due 
to pions from resonance decays (e.g. p, A), although such 
an effect is reduced by the lower transverse momentum 
cut. An alternative explanation is that hydrodynamical 
collective behavior may develop even in the small p + p 
system, which we will explore at Section [V Al 

In Fig. [TUl the collision energy dependence of Tiny is 
shown by compiling results from past experiments |24[ . 
m, [48l - [5g | . The values of Tjnv reported here are obtained 
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FIG. 9: Inverse slope parameter Tinv for tt*, , p and p in 
p + p collisions at ^ = 200 GeV and 62.4 GeV. The fitting 
range is ttit — m — 0.3-1.0 GeV/c^ for all particle species at 
both collision energies. The errors are statistical and system- 
atic combined. The statistical errors are negligible. 



TABLE V: Fitting results for A, Tinv of Eq.fTOlfor n^, K"^, p, 
and p in 200 and 62.4 GeV p + p collisions. The fitting range 
is fixed as ttit — m = 0.3-1.0 GeV/c^ for all particle species 
at both collision energies. 





hadron 


A 


21iiv 


xVndf 


(GeV) 






(GeV/c2) 




200 


7r+ 


73.4 ±7.1 


0.190 ± 0.005 


5.6/5 




7r~ 


74.8 ± 7.2 


0.189 ± 0.005 


3.1/5 




K+ 


3.25 ± 0.29 


0.232 ± 0.007 


3.6/6 




K- 


2.99 ± 0.27 


0.239 ± 0.008 


3.6/6 




P 


0.85 ± 0.14 


0.245 ± 0.014 


1.0/7 




P 


0.74 ± 0.13 


0.241 ± 0.014 


0.5/7 


62.4 


7r+ 


61.7 ± 5.9 


0.182 ± 0.005 


3.1/5 




7r~ 


65.2 ± 5.3 


0.179 ± 0.004 


4.7/5 




K+ 


2.44 ± 0.22 


0.219 ± 0.007 


2.6/6 




K- 


2.21 ± 0.20 


0.213 ± 0.006 


4.6/6 




P 


0.81 ± 0.10 


0.227 ± 0.010 


1.1/7 




P 


0.49 ± 0.07 


0.221 ± 0.010 


0.3/7 



by fitting all the pr spectra in the same way. The fitting 
range is mT — m — 0.3-1.0 GeV/c^ for all particle species 
in all collision systems. The T-^v values for RHIC ener- 
gies are consistent with earlier experimental results at 
other energies [24, 48-50]. For both pious and kaous, the 
inverse slope parameters increase with collision energy 
from Tinv = 120 MeV/c^ to 170 MeV/c^ (240 MeV/c^) 
for pions (kaons) up to y^s = 200 GeV. According to 
Tevatron data, Tinv seems to be saturated at i/s above 
200 GeV. The inverse slope parameters of protons and 
antiprotons indicate an increase at lower ^/s which keeps 
on increasing even at Tevatron energies. We look forward 
to data from the LHC to further clarify these issues. 



D. Particle ratios 



The 'K~ /-K^ and K~ jK^ ratios show a flat pt depen- 
dence at both 200 and 62.4 GeV energies. The 7r~/7r"'" 
ratio is almost unity at both energies. The / ratio 
is consistent with unity at ^/s = 200 GeV, while it de- 
creases to 0.8-0.9 in the measured pt range at 62.4 GeV. 
On the other hand, p/p ratio seems to be a decreasing 
function of pT at 200 GeV, from the value of ~ 0.8 at 
PT = 1.0 GeV/c to 0.6 at pt = 4.5 GeV/c. Note that 
we fit the p/p ratio for 200 GeV p + p from px = 1-4.5 
GeV/c to a linear function, a -\- bpT, which gives a = 
0.93 ± 0.02, b = -0.07 ± 0.01. This decrease, also seen 



at lower ^/s [2J], might be the result of a difference of 
fragmentation between quark jet and gluon jet at high 
Pt re gion as suggested by the DSS fragmentation func- 
tions [47'|. However, the NLO pQCD calculation using 
the DSS fragmentation functions (lines on the panels for 
p/p ratios) shows this effect is in disagreement with the 
measured p/p ratios. At 62.4 GeV, we can not conclude 
on the significance of the decrease of the p/p ratios as a 
function of pt due to large statistical fluctuations. It is 
important to note agreement of the ISR measurements of 
the anti-particle to particle ratios as a function of pt at 
y/s = 62 A GeV (Fig. [TT|) with the present measurements 
except for the p/p ratio where there is a large discrepancy. 
The p/p ratio integrated over all pt decreases from 0.8 
at 200 GeV to 0.5 at 62.4 GeV (see further discussion 
in Section HVEp . At low pT, the large systematic uncer- 
tainties of the p/p ratio are due to the uncertainties of 
the weak decay feed down corrections. 

Figure [12] presents the ratios of K^/tt^, K^/tt^, 
p/7r+, p/tt", p/tt", and p/n° as a function of pt- Both 
the /tt^ and the K~ /tt~ ratios increase with increas- 
ing Pt up to the pt = 2 GeV/c limit of the measurement. 
Both the p/tt" and the p/tt^ ratios seem to increase with 
PT for PT > 2 GeV/c, although the p/7r° ratio is rela- 
tively flat at ^/s — 200 GeV in the same transverse mo- 
mentum region. Clearly, better statistics are required to 
reach a firm conclusion. As a function of y/s the /tt^, 
p/tt^ and p/ir^ ratios do not change significantly, while 
the K~ /Tr~ ratio increases and the p/tt"*" andp/Tr*^ ratios 
decrease significantly for pT > 1 GeV/c as the collision 
energy is increased from = 62.4 GeV to 200 GeV. 



E. 



r) and dN/dy 



Figures [TT] and [T2] show particle ratios such as 
antiparticle-to-particle, K/tt, p/tt as a function of px- 



Mean transverse momentum (py) and particle yield 
per unit rapidity dN/dy are determined by integrating 
the measured px spectrum for each particle species. For 
the unmeasured pT region, we fit the measured pt spec- 
trum with a Tsallis function (52l | given below, as in a 
related publication j28| . and also with an ttit exponen- 
tial function, then extrapolate the obtained function to 
the unmeasured pr region. The pr ranges for fitting are: 
0.4-3.0 GeV/c for pions, 0.4-2.0 GeV/c for kaons, 0.5-4.0 
GeV/c for protons and antiprotons. 

The final yield dN/ dy is calculated by taking the sum 
of the yield from the data, and the yield from the func- 
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FIG. 10: (color online) Inverse slope parameter Tinv for n'^ + n , K'^ + K , p and p in p + p and p + p collisions versus 
collision energy ^/s. The fitting range is niT — m = 0.3-1.0 GeV/c^ for all particle species for all collision systems. The 
errors are statistical and systematic combined. The statistical errors are negligible. The other experimental data are taken 
from 113, [3^, ^8,-^]. 



TABLE VI: Fitting 


results from using 


the Tsallis distribution (Eq.[TTJ for tt*, 


p and p in p - 


-p collisions at ^/s 


= 200 and 


62.4 GeV. 














(GeV) 


hadron 


dN/dy 


q 




C 


xVndf 


200 


7V+ 


0.963 ± 0.071 


8.24 ± 0.33 


0.115 


± 0.006 


4.3/23 




7r~ 


0.900 ± 0.063 


8.95 ± 0.39 


0.123 


± 0.006 


3.2/23 




K+ 


0.108 ± 0.006 


6.25 ± 0.64 


0.137 


± 0.011 


1.6/13 




K- 


0.103 ± 0.005 


7.00 ± 0.78 


0.147 


± 0.011 


2.9/13 




P 


0.044 ± 0.004 


11.1 ± 1.6 


0.184 


± 0.014 


4.1/22 




P 


0.037 ± 0.003 


12.0 ± 1.8 


0.186 


± 0.014 


1.3/22 


62.4 


n+ 


0.782 ± 0.056 


12.1 ± 0.9 


0.133 


± 0.007 


4.6/22 




7r~ 


0.824 ± 0.053 


11.9 ± 0.7 


0.128 


± 0.006 


4.8/22 




K+ 


0.076 ± 0.003 


10.2 ± 1.8 


0.165 


± 0.012 


4.9/13 




R- 


0.067 ± 0.003 


11.6 ± 2.1 


0.164 


± 0.011 


2.2/13 




P 


0.040 ± 0.003 


24.5 ± 9.9 


0.201 


± 0.015 


7.1/21 




P 


0.022 ± 0.002 


32.5 ± 21.0 


0.202 


± 0.018 


7.9/21 



tional form in the unmeasured px region. The total in- 
elastic cross sections are assumed to be 42.0 mb and 35.6 
mb for 200 GeV and 62.4 GeV respectively. For (pt), we 
integrate the measured px spectrum with px weighting, 
and then divide it by the obtained dN/dy. The final val- 
ues are obtained by averaging the results of the two fits. 
The systematic uncertainties are evaluated as half of the 
difference between these fitting values. 

• Tsallis distribution is given by Eq. [TT] below. In 
this fitting form, the free parameters are dN/dy, 
q and C, while the mass m is fixed to the hadron 
mass. The fitting results are given in Table IVll 



1 



d^N 



2iipT dydpT 



dN 



(Q-l)(<7-2) 



dy 2nqC[qC + m{q-2)] 



rriT — m 



qC 



(11) 



Exponential distribution in tot is given by Eq. [T^] 
below. The free fit parameters are the normaliza- 
tion constant A and the inverse slope parameter 



1 d^N 
2ttpt dydpT 



Aexp(-^) 

^ ^ inv ^ 



dy 



(12) 



(13) 
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FIG. 11: (color online) -k~ /n'^ , K~ /K^ , p/p ratios as a func- 
tion of pt in p+p collisions at -^s — (left) 200 and (right) 62.4 
GeV. Systematic uncertainties are shown as vertical shaded 
band. STAR data are from [5l[ and ISR data are from [2i |. 
For p/p ratios, the NLO pQGD calculations using the DSS 
fragmentation functions fiTil are also shown as (solid lines) 
fj. = Pt and (dashed lines) fi — 2pT, Pt/2. 
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The obtained (px) values are summarized in Table IVlIl 
They are plotted in Fig. [T31 which indicates a clear in- 
crease of {pt) with hadron mass. The values at 200 GeV 
are almost same as those for 62.4 GeV data. If the spec- 
tral shape is a pure exponential, (pt) should be equal to 
2Tinv analytically. By comparing Tables |V] and I VIII the 
measured {pr) is almost 2Tinv for pions. But for kaons 
and (anti)protons, the measured (pt) is systematically 
larger than 2 Tinv This demonstrates that the spectral 
shape at low px is not a pure exponential especially for 
kaons and (anti)protons. 

The collision energy dependence of (jjt) for each parti- 
cle type is shown in Fig. [TH Data shown here are: lower 
energy data [H], ISR data Tevatron data [H, HI], 
and RHIG data from STAR and PHENIX (present 
study). The (px) values for all the other experiments 
have been determined by fitting the px spectra. For pi- 



FIG. 12: (color online) (top) K~^/tv~^ and K~ /it~ ratios, 
(middle) p/tt"*" and (bottom) p/vr^ ratios, and (bottom) p/n^'' 
and P/tt" ratios as a function of pt in p-l-p collisions at ^/s — 
(left) 200 and (right) 62.4 GeV. Systematic uncertainties are 
shown as vertical shaded bands. 



ons, the (px) shows a linear increase in ln{^/s). For kaons 
and (anti)protons the increase is much faster than those 
for pions. However, systematic issues at both lower and 
higher center-of-mass energies remain to be resolved. 

Figure [TSl shows the dependence of (pr) on the central- 
ity of the collisions (given by the number of participating 
nucleons, A^part) for tt^, K^, p and p in Au-|-Au coUisions 
at y/SNN = 200 GeV [33| as compared to minimum bias 
p + p collisions at ^/s = 200 GeV (present analysis) . The 
error bars in the figure represent the statistical errors. 
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The systematic errors from cut conditions are shown as 
shaded boxes on the right for each particle species. The 
systematic errors from extrapolations, which are scaled 
by a factor of two for clarity, are shown in the bottom for 
each particle species. It is found that {pt) for all particle 
species increases from the most peripheral to midcentral 
collisions, and appears to saturate from the midcentral 
to central collisions. The (pt) inp+p are consistent with 
the expectation from the A^part dependence in Au+Au, 
and are similar to the values in peripheral Au+Au. 
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FIG. 13: Mean transverse momentum (pr) as a function of 
mass in p+p collisions at y^s = 200 and 62.4 GeV. The errors 
are statistical and systematic combined. The statistical errors 
are negligible. 



TABLE VII: Mean transverse momentum {{pt)) and parti- 
cle yield (dN/dy) for tt"^ , , p and p in p + p collisions 
at — 200 and 62.4 GeV. The errors are statistical and 
systematic combined, but the statistical errors are negligible. 





hadron 


(pt) 


dN/dy 


(GeV) 




(GeV/c) 





n+ 


0.379 


± 





021 


0.842 


± 


0.127 




0.385 


± 





014 


0.810 


± 


0.096 


K+ 


0.570 


± 





012 


0.099 


± 


0.010 


R- 


0.573 


± 





014 


0.096 


± 


0.009 


P 


0.696 


± 





025 


0.043 


± 


0.003 


P 


0.698 


± 





023 


0.035 


± 


0.002 


7r+ 


0.373 


± 





013 


0.722 


± 


0.066 




0.366 


± 





016 


0.750 


± 


0.079 


K+ 


0.558 


± 





012 


0.072 


± 


0.004 


K- 


0.544 


± 





013 


0.064 


± 


0.004 


P 


0.710 


± 





023 


0.034 


± 


0.002 


P 


0.709 


± 





040 


0.018 


± 


0.001 



Figure [17] shows dN/dy as a function of collision energy 
for each particle species. Our results on dN/dy are con- 
sistent with the those at ISR energies It should be 
noted that STAR quotes the nonsingle diffractive (NSD) 
multiplicity while our measurement quotes the inelastic 
multiplicity, normalizing the integrated measured inclu- 
sive cross section by the total inelastic cross section fs^ . 
At ^/s = 200 GeV, the inelastic cross section [a^^^^) is 
42 mb [5^, and the single diffractive (SD) cross section is 
almost equal to the double diffractive (DD) cross section, 
^fm ~ ^NN ~ 4 mb [11]. As the single diffractive refers 
only to the projectile proton in a p -f p fixed target mea- 
surement, one has to subtract the SD cross section for 
each proton from the inelastic cross section to determine 
the NSD cross section [s^]. The resulting NSD cross sec- 
tion (cr^-5-°) should be 42 - 2 X 4mb = 34 mb. The ratio 
of the NSD multiplicity to the inelastic multiplicity is: 
^INEL/^NSD ^ 42/34 = 1 24, i.e. the NSD multiplic- 
ity is 24% higher than the inelastic multiplicity, and this 
effect can be actually seen in the experimental data [5^ . 

We would like to point out also that the NSD charged 
particle muhiphcity at = 200 GcV by STAR is « 
20% larger than other NSD resuhs [58]. By taking this 
fact and the difference between NSD and inelastic into 
account, one can naturally understand « 50% difference 
on yields between STAR and the present analysis, for 
pious and kaons, as shown in Fig. 1171 For protons and 
antiprotons the difference between STAR and the present 
analysis is bigger than those in pious and kaons. In ad- 
dition to the effects we have mentioned above, the weak 
decay feed-down correction can contribute to it, since we 
remove p and p from the weak decay (see Section IlII Dp . 
while STAR does not. 

Figure [18] shows the collision centrality dependence of 
dN/dy per participant pair (0.5 iVpart) in p + p (present 
analysis) and Au-|-Au [33 collisions at ^/sWn — 200 GeV. 
The error bars on each point represent the quadratic sum 
of the statistical errors and systematic errors from cut 
conditions. The statistical errors are negligible. The 
lines represent the effect of the systematic error on A^part 
which affects all curves in the same way. The data in- 
dicate that dN/dy per participant pair increases for all 
particle species with A^part up to ~ 100, and saturates 
from the midcentral to the most central collisions. As 
seen in Fig. [15] for (pr), the dN/dy values in p + p are 
consistent with the expectation from the iVpart depen- 
dence in Au-|-Au. 



The dN/dy values at midrapidity are summarized in 
Table IVIII They are plotted in Fig. [16] as a function of 
hadron mass for both 200 and 62.4 GeV collision energies. 
There are differences in the yield between 200 and 62.4 
GeV especially for kaons and antiprotons continuing the 
trend observed at lower y/s [241]. It is interesting to note 
that even in the situation that dN/dy is different between 
^/s = 200 and 62.4 GeV, {pt) is quite similar for both 
energies. 



F. Nuclear Modification Factor Raa 

In order to quantify the modification effect in nucleus- 
nucleus (A+A) collisions with respect to nucleon-nucleon 
collisions, the nuclear modification factor i?AA is used. 
i?AA is the ratio between the yield in A-l-A scaled by 
the average number of binary nucleon-nucleon collisions 
((A"cou)) and the yield in p+p, as defined by the following 
equation: 
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FIG. 14: (color online) Mean transverse momentum (pt) for 7r+ + ir , K'^ + K , p and p as a function of ^/s in p + p and 
p + p collisions [2^. [3^. IZsl. l49l. Issl] . The errors are statistical and systematic combined. The statistical errors are negligible. 
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FIG. 16: Particle yield dN/dy as a function of mass in p + p 
collisions at y/s = 200 and 62.4 GeV. The errors are statistical 
and systematic combined. The statistical errors are negligible. 
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FIG. 15: (color online) Mean transverse momentum as a 
function of centrality (A'^part) for pious, kaons, protons and 
anti-protons in in p + p (present analysis, red color) and 
Au+Au [131 (black color) collisions at ^/sWh = 200 GeV. 
The left (right) panel shows the (pt) for positive (negative) 
particles. The error bars are statistical errors. The system- 
atic errors from cuts conditions are shown as shaded boxes 
on the right for each particle species. The systematic errors 
from extrapolations, which are scaled by a factor of two for 
clarity, are shown in the bottom for protons and anti-protons 
(dashed-dot lines), kaons (dotted lines), and pions (dashed 
lines) . 



Raa{pt) 



(l/jVrA) d^NAA/dpTdy 
(Taa) X dL^app/dprdy 



(14) 



where (Taa) is the nuclear thickness function, defined as 
follows: (Taa) = (iVcoii)/crj,p°'. For the total A+A inter- 
action cross section ct^^ (minimum bias A-l-A collisions). 



(rAA)-AV4l. ^ 

In general, i?AA is expressed as a function of pT and 
collision centrality for A-l-A collisions. Due to the dom- 
inance of hard scatterings of partons at high pT, Raa 
is expected to be around unity above pt ~ 2 GeV/c, if 
there is no yield modification by the nucleus in A-l-A. If 
there is a suppression (enhancement), i?AA is less than 
(greater than) unity. For the total A-l-A interaction cross 
section at a given pT integrated over centrality (mini- 
mum bias A-l-A collisions) ctaa(pt) — A^o-ppipr) and 
Raa{pt) = 1.0. 

Figure [12] shows the i?AA of n^, tt" , K"^ , p, and p in 
Au+Au collisions at y^SNN = 200 GeV at 0-5 % colli- 
sion centrality. The data for identified charged hadrons 
in Au+Au at ^/snn — 200 GeV are taken from mea- 
sured by the PHENIX experiment, and those for p+p are 
taken from the present analysis at ^/s = 200 GeV. The 
i?AA for neutral pions is taken from [s^. The overall 
normalization uncertainty on i?AA (13.8%) is shown as a 
shaded box around unity (at = 0.1 GeV/c)), which 
is the quadratic sum of (1) the uncertainty oi p + p in- 
elastic cross section (9.7%) and (2) the uncertanty (A^coii) 
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FIG. 17: (color online) (upper) Particle yield dN/dy at midrapidity for (tt^ +-k )/2, (K^ +K )/2 as a function of in p + p 
collisions [23, iS]. The errors are statistical and systematic combined, but the statistical errors are negligible. The dN/dy from 
STAR is determined for NSD p + p events, (lower) Similar plots for p and p with feed-down correction applied to our data. 
The dN/dy from STAR is determined for NSD p + p events, and is not corrected for weak-decay feed down. 



(9.9%). 



V. DISCUSSION 



For pions i?AA is largely suppressed by a factor of ss 5, 
compared to p + p. This suppression effect is understood 
to be due to jet quenching or energy loss of partons in hot 
and dense medium created in Au-|-Au central collisions at 
RHIC energies [13, HH . For kaons there is a similar trend 
as pions over a more limited px range. For protons and 
antiprotons there is an enhancement in pT — 2-4 GeV/c. 
As reported in [2^ [s^, Isl] , possible explanations of the 
observed enhancements include the quark recombination 
model [63l-l65l| and/or strong partonic and hadronic radial 
flow tef 



In this section, we discuss (1) soft particle production 
at low pt, including the possibility of radial flow in p + p 
collisions, and (2) the transition from soft to hard pro- 
cess, and hadron fragmentation at high pT, where we 
show the XT scaling of measured spectra, and make a 
comparison with NLO and NLL pQCD calculations. 



A. Radial flow 

In heavy ion collisions at RHIC energies, it is found 
that the inverse slope parameter (Tinv) of rriT — m spectra 
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FIG. 18: (color online) Particle yield per unit rapidity 
(dN/dy) per participant pair (0.5 A'^part) as a function of 
A'^part for pions, kaons, protons and antiprotons in p + p 
(present analysis, red color) and Au+Au collisions [37[ at 
yfs]^ = 200 GeV. The left (right) panel shows the dN/dy 
for positive (negative) particles. The error bars represent the 
quadratic sum of statistical errors and systematic errors from 
cut conditions. The lines represent the effect of the systematic 
error on A'part which affects all curves in the same way. 
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FIG. 19: (color online) Raa of tt^, tv" , K^, p, and p in 
Au+Au collisions at y'sjviv = 200 GeV at 0-5 % collision 
centrality. The data for identified c harg ed hadrons in Au+Au 
at y^sjvjv = 200 GeV is taken from ^3 a-nd the data for p+p 
is taken form the present analysis at y/s = 200 GeV. The 
neutral pions (PHENIX) is taken from [5^. The statistical 
uncertainties are shown as bars, and the systematic uncer- 
tainties are shown as shaded boxes on each data point. The 
overall normalization uncertainty on -Raa (13.8%) is shown 
in the shaded box around unity (at pt = 0.1 GeV/c)), which 
is the quadratic sum of (1) uncertainty of p + p inelastic cross 
section (9.7%) and (2) uncertainty (iVcoii) (9.9%). 



has a clear dependence on the hadron mass, i.e. heavier 
particles have larger inverse slope parameters (stI . [67j . 
Tiiiv increases almost linearly as a function of particle 
mass; is largest vifhen the nucleus-nucleus collision has a 
small impact parameter (central collisions); and is small- 
est for the collisions with a large impact parameter (pe- 
ripheral collisions), as shown in Fig. 

This experimental observation can be interpreted as 
the existence of a radial flow generated by violent 
nucleon-nucleon collisions in two colliding nuclei and de- 
veloped both in the quark-gluon plasma (QGP) phase 
and in hadronic rescatterings [i^ . The radial flow veloc- 
ity increases the transverse momentum of particles pro- 
portional to their mass, thus Tinv increases linearly as a 
function of particle mass. It is interesting to determine 
whether or not such an expansion is observed in high 
energy p + p collisions 53] . 

Figure [20] shows the mass dependence of inverse slope 
parameter Tjnv in rriT — m spectra for positive (left) and 
negative (right) particles in p + p collisions at ^/s = 200 
and 62.4 GeV (also shown in Fig. IH), as well as for pe- 
ripheral, midcentral and central in Au+Au collisions at 
./sJTn = 200 GeV [STj. The fit ranges are tot-w = 0.2- 
1.0 GeV/c^ for pions, and mr — m = 0.1-1.0 GeV/c^ for 
kaons, protons, and antiprotons, which are chosen in or- 
der to perform a fair comparison with Tjnv in Au+Au 
collisions at RHIC (sj • The values of Tinv in p + p for 
these fit ranges (see Table IVIIip are all lower by roughly 
one standard deviation from the values in Table |V] for 
the common fit range of tot — m = 0.3-1.0 GeV/c^. 
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FIG. 20: (color online) Mass dependence of inverse slope pa- 
rameter Tinv in mr — w- spectra for (left) positive and (right) 
negative hadrons in p+p collisions at y^s = 200 and 62.4 GeV, 
as well as for peripheral, midcentral and central in Au+Au 
collisions at ^Jsnn = 200 GeV [37[. The errors are statistical 
and systematic combined, smaller than symbols. The statis- 
tical errors are negligible. The fit ranges are 0.2-1.0 GeV/c^ 
for pions and 0.1-1.0 GeV/c^ for kaons, protons, and antipro- 
tons in rriT — ra [s^. The dotted lines represent a linear fit 
of the results for each data set, as a function of mass using 
Eq. 

In general, the inverse slope parameters increase with 
increasing particle mass in both Au+Au and p + p colli- 
sions at 200 GeV. However, this increase is only modest in 
p+p collisions and slightly weaker than in 60-92% central 
Au+Au collisions at 200 GeV. Also note that there is a 
mean multiplicity dependence of the transverse momen- 
tum spectra in p + p collisions [1^ , that is not discussed 
in the present paper. 
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TABLE VIII: Inverse slope parameter Tinv for tt*, , p and 
p for p + p collisions at = 200 and 62.4 GeV. The fit 
ranges are 0.2-1.0 GeV/c^ for pions and 0.1-1.0 GeV/c^ for 
kaons, protons, and antiprotons in mr — m. These fit ranges 
are chosen in order to perform a comparison with Tinv in 
Au+Au collisions at RHIC [S^]. The errors are statistical and 
systematic combined, but the statistical errors are negligible. 
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IX: The extracted fit parameters of the freeze-out 


temperature (To) and the measure of the strength of the av- 


erage radial transverse flow {{ut)) using Eq. (|15|) 
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We use a radial flow picture [H, with the fitting 
function: 

T = ^o^-m(^it)^ (15) 

where To is a hadron freeze-out temperature and {ut) 
is a "measure" of the strength of the (average radial) 
transverse flow. The relationship between the averaged 
transverse velocity ((/3t)) and {ut) is given by 

(A) = {m)/^jl + {utf. (16) 

The dotted lines in Fig. [20] represent the linear fit to 
the p + p collisions at -^s = 200 and 62.4 GeV, that are 
compared to those in Au-|-Au collisions at y/s^ = 200 
GeV in three different collision centrality classes. The fit 
results in p + p are also given in Table IIXI For Au-|-Au 
most central data (0-5%), (ut) w 0.49 ±0.07, while in 
p + p, {ut) « 0.28 at both 62.4 GeV and 200 GeV. While 



this radial flow model is consistent with the data in cen- 
tral and midcentral Au-|-Au, i.e. n/K/p points are on 
a straight line, it does not give a good description of ei- 
ther peripheral Au-f Au or p + p collisions (poor i^i 
Table IIX[) . This can be interpreted that radial fiow is 
absent in p + p, where the tt/ K/p points are obviously 
not on a straight line (Fig. [S]), and that the radial flow 
develops only for a larger system. 

B. XT scaling 

From the measurements of px spectra of hadrons in 
p + p collisions, it is known that fragmentation of hard 
scattered partons is the dominant production mecha- 
nism of high pt hadrons. It has been predicted the- 
oretically from general principles, that such a produc- 
tion mechanism leads to a data scaling behavior called 
"xt scaling" , where the scaling variable is defined as 
XT — "ipT I y/s- Such a data scaling behavior was seen first 
on preliminary ISR data at CERN as reported in [36| . 

In the kinematic range corresponding to the xt scaling 
limit, the invariant cross section near midrapidity can be 
written as: 

where F{xt) and G{xt) are universal scaling functions. 
The parameter ricff is characteristic for the type of in- 
teraction between constituent partons. For exarnple for 
single photon or vector gluon exchange, = 4 [l|. Be- 
cause of higher order effects, the running of the strong 
coupling constant as = as{Q^), the evolution of the par- 
ton distribution functions and fragmentation functions, 
and nonvanishing transverse momentum kx of the initial- 
state, ricff in general is not a constant but a function of xt 
and y/s, i.e. rioff = rioff (a;T, v^)- This rics corresponds 
to the logarithmic variation of yield ratios at the same 
Xt for different y/s [t^. Note that the xt scaling power 
rice is different from the exponent n that characterizes 
the power-law behavior of the single particle invariant 
spectrum at high pT ■ 

The value of rics depends on both the value of ^/s and 
the range of xt and depending on the reaction tends to 
settle at an asymptotic value between 6 and 4.5 where 
hard-scattering dominates and higher twist effects are 
small. This fact can also be used to determine the tran- 
sition between soft and hard particle production mecha- 
nisms. 

Earlier measurements of tLcs (xt , in P + P collisions 
found values in the range of 5-8 [3^, [ZlnZJI • Here we 
present the PHENIX results for the xt scaling of pions, 
protons, and antiprotons and compare them with earlier 
data measured at various different values of ^/s. Due to 
the limited pt range of our kaon measurements, kaons 
are not included in these comparisons. 

We have evaluated the XT-scaling power rioff using two 
different methods that are both based on Eq. ([T7|) : 
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FIG. 21: (color online) (upper left) a:T-scaling power riofr as determined from the ratios of yields as a function of xt, for (open 
circle) neutral pions, (open square) protons, and (open triangle) antiprotons using p + p data at -^s — 200 and 62.4 GeV 
energies. The error of each data point is from the systematic and statistical errors of pr spectra. The other plots show xt 
spectra for (lower left) pions (tt^, tt"), (upper right) protons, and (lower right) antiprotons in p + p collisions at different y^s 
at midrapidity. Only statistical uncertainties are shown. The dashed curves are the fitting results. 



Method 1 is based on the hnear variation of the loga- 
rithm of the ratio of the yields at different y/s: 



\og{Yield(xT, 62 A) /Yield{xT, 200)) 
log(200/62.4) ' 



(18) 



The ncs{xT) is shown in Fig. [21] as a function of 
Xt for neutral pions, protons, and antiprotons for 
p + p coUisions at RHIC. 

Method 2 is based on fitting the xt distributions for a 
given type of particle measured at different ener- 
gies. A common fitting function is defined as fol- 
lows: 



E- 



XT 



(19) 



limiting the fitting region to the high transverse 
momentum region (jjt > 2 GeV/c). 

The Xt distributions for pions, protons, and antipro- 
tons are shown in Fig. [5T] PHENIX data are presented 



together with earlier data of (2J, |30|, |3J, l75|. Dashed 
curves show the fitting results. The obtained rics values 
are summarized in Table [Xj 

The exponent riofi of the xt scaling is found to have 
similar values for different particles, in the range of 6.3- 
6.5 for pions, protons, and antiprotons. The data points 
deviate from the xt scaling in the transverse momentum 
region of < 2 GeV/c. This scaling violation may be 
interpreted as a transition from hard to soft multipar- 
ticle production. For the highest xt points for protons 
and antiprotons (but not for pions) the asymptotic xt 
curve gets steeper. Further measurements at larger xt, 
possibly at lower center-of-mass energies are needed to 
clarify this point. 



Comparison to NLO and NLL pQCD 
calculations 



In Figs. [22] to nil our results for pion, proton and an- 
tiproton pT spectra at ^/s = 200 and 62.4 GeV in p + p 
collisions are compared to the NLO pQCD calculations. 
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FIG. 22: (color online) Transverse momentum distributions 
for (upper) positive and (lower) negative particles at ^ — 
200 GeV in p + p collisions. Only statistical uncertainties are 
shown. The normalization uncertainty (9.7%) is not included. 
NLO pQCD calculations (FF: DSS) are also shown. Solid 
lines are for fi = pT, and dashed lines are for fi = pt/2, 2pT. 
The lower panel in each plot shows the ratio of (data - pQCD 
result)/pQCD result for each particle species. 



TABLE X: Summary of a;T-scaling power Ticft in p + p colli- 
sions. The errors are systematic error from the fitting. 
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xVndf 
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0.84 ± 0.04 
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Because of the limited px reach in the measurements, the 
results for charged kaons are not compared to the NLO 
pQCD calculations. In these NLO pQCD calculations 
for 77 < 1 from W. Vogelsang [43| the cross section is fac- 
torized into initial parton distribution functions (PDFs) 
in the colliding protons, short-distance partonic hard- 
scattering cross sections which can be evaluated using 
perturbative QCD, and parton-to-hadron fragmentation 
functions (FFs). 



For the description of the initial parton distributions, 
the Coordinated Theoretical-Experimental Project on 
QCD (CTEQ6M5) [76] PDFs are used. Different scales 
/i = Pt/2, pt, '^Pt are utilized, where fj. represents equal 
factorization, renormalization, and fragmentation scales. 
These provide initial conditions for the pQCD cross 
section calculations. Partons are then fragmented to 
hadrons with the help of the de Florian-Sassot-Stratmann 
(DSS) set of fra gme ntation functions which have the 
charge separation fTTi]. There are several other FFs, such 
as the "Albino-Knichl-Kramer" (AKK) set of FFs W^, 
and the "Kniehl-Kramer-Potter" (KKP) set of FF UM- 



Only the results for DSS FFs are shown in this paper, 
because they give better agreement with our measure- 
ments than other FFs. For example, in p + p collisions at 
= 200 GeV the yields for {p + p)/2 in AKK (KKP) 
FFs are a factor of two smaller (larger) than the present 
measurement. 

It is known that pion production in ^/s = 200 GeV p+p 
collisions is reasonably well described by pQCD down to 
~ 2 GeV/c and up to pr ~ 20 GeV/c P, On 
the other hand, there are large variations in the p and p 
yields among various fragmentation functions |33| as we 
mentioned above. From the comparisons between baryon 
data and pQCD calculations at both ^/s = 200 and 62.4 
GeV, it is potentially interesting to obtain a constraint on 
the fragmentation function, particularly the gluon frag- 
mentation function for p and p. 

For the DSS fragmentation function, there is good 
agreement between the data and NLO pQCD calcula- 
tions for pious and protons at 200 GeV, but not so good 
agreement with p. It is more clearly shown in Fig. 1111 
that the p/p ratio at 200 GeV is not correctly described 
with NLO + DSS framework, which indicates that there 
is still room to improve the DSS fragmentation functions. 
The left-side plots of Fig. US] show that for 62.4 GeV NLO 
+ DSS pQCD calculations underestimate yields by a fac- 
tor of two or three for all species. However, as it is still 
on the edge of the scale uncertainty of the NLO calcula- 
tion, NLO pQCD agrees with the data within the large 
uncertainties. 

As shown in [34| , the NLL calculations give much bet- 
ter agreement with the data for tt" in p + p collisions at 
•y/s = 62.4 GeV. This means the resummed calculation 
is necessary to describe the cross section at 62.4 GeV. 
On the other hand, the resummation for y/s = 200 GeV 
is not reliable, since the resummation can be done for a 
larger xt — 2pT / which is not accessible for ^/s — 200 
GeV data due to the pr limitation. The right-side plots 
of Fig. [23] show the px distributions for tt*, p, and p in 
p + p collisions at 62.4 GeV, together with the results of 
NLL pQCD calculations. The DSS FFs are used. It is 
found that the agreement between NLL pQCD and data 
is better than those for NLO pQCD. 

The presented px spectra extend to the semi-hard 3-4 
GeV/c region for pions and (anti)protons, which make 
them useful as a baseline to study in further detail the 
nuclear modification factor in A-l-A collisions. More de- 
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FIG. 23: (color online) Transverse momentum distributions for (upper) positive and (lower) negative particles at ^/s = 62.4 
GeV in p + p collisions. Only statistical uncertainties are shown. The normalization uncertainty (11%) is not included, (left) 
NLO and (right) NLL pQCD calculations (FF: DSS) are also shown. Solid lines are for /j, = pr, and dashed lines are for /i — 
Pt/^, 2pT- The lower panel in each plot shows the ratio of (data - pQCD result)/pQCD result for each particle species. 



tailed measurements at larger pT are necessary for the 
further understanding of FFs and its particle species de- 
pendence at each beam energy. 

VI. SUMMARY AND CONCLUSION 

We presented transverse momentum distributions and 
yields for ir^, K^, p and pin p+p collisions at ^/s = 200 
and 62.4 GeV at midrapidity, which provide an important 
baseline for heavy- ion-collision measurements at RHIC. 
The inverse slope parameter Tjnv, mean transverse mo- 
mentum (pr) , and yield per unit rapidity dN/dy are com- 
pared to the measurements at different ^/s in p -\- p and 
p + p collisions. While Tinv and {pr) show a similar value 
for all particle species between 200 and 62.4 GeV, dN/dy 
shows a relatively large difference, especially for kaons 
and antiprotons, between 200 and 62.4 GeV. The p/p ra- 
tio is -0.8 at 200 GeV, and -0.5 at 62.4 GeV and the 
Pt dependence of the p/ir^ {p/ir'^) ratios varies between 
62.4 and 200 GeV. Together with the measured dN/dy 
this gives insight on baryon transport and production at 
midrapidity. 



We also analyzed the scaling properties of identified 
particle spectra, such as the mx scaling and xt scaling. 
Baryons and mesons are split in the mx spectral shape 
at both 200 and 62.4 GeV. This splitting can be under- 
stood as the difference of hard production yields between 
baryons and mesons. The xt scaling power rioff shows 
similar values for pions, protons, and antiprotons. 

We also compared the results in p + p collisions at 
200 GeV with the those in Au-hAu at 200 GeV in the 
same experiment. It is found that Tjnv, {pt), and dN/dy 
change smoothly from p + p to Au-t-Au, and all the val- 
ues in p + p are consistent with expectations from the 
-^part dependence in Au+Au. On the nuclear modifica- 
tion factor i?AA, there is a large suppression for pions, 
while there is an enhancement for protons and antipro- 
tons at Pt =2-4 GeV/c. The observed suppression can 
be understood by the energy loss of partons in hot and 
dense medium created in Au-|-Au central collisions at 
RHIC energies [g^, HH . Possible explanations of the ob- 
served enhancements for pro tons and antiprotons include 
quark recombination l63l - l65| and/or strong partonic and 
hadronic radial flow 66] . 

Identified particle spectra are extended to the semi- 
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hard 3 4 GcV/c region for pions and (anti)protons, which 
make it possible to study in further detail the nuclear 
modification factor of identified particles in A+A colli- 
sions. NLO pQCD calculations with DSS fragmentation 
function show good agreement for pions and protons at 
200 GeV while there is less good agreement with p. This 
indicates that fragmentation functions can be improved 
further. 

For 62.4 GeV, NLO pQCD calculations underestimate 
by a factor of two or three the yields for all particle 
species. On the other hand, NLL pQCD calculation gives 
a better agreement with the data. This suggests that re- 
summed calculations are necessary to describe the cross 
section at 62.4 GeV. 

From comparisons to some calculations such as 
NLO/NLL pQCD framework, one can discuss the mech- 
anism of soft and hard particle production m p + p col- 
lisions. There is a transition between these two regions 
("soft-hard transition") at pr ~2 GeV/c for pions, and at 
Pt ^3 GeV/c for (anti)protons, or equivalently, m,T — m, 
= 1-2 GeV/c^ for all particle species at both energies. 
The fractions of soft and hard components gradually 
change in the transition region. The new measurements 
presented in this work indicate that the behavior of T^nv 
and {pt) of identified particles in p -|- p collisions needs 
improved systematics of measurement at both higher and 
lower ^/s to clarify the precise \fs dependences. 
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Appendix A: Table of cross sections 



The cross sections for tt^ , , p and p in p+p collisions 
at -^/s = 200 and 62.4 GeV at midrapidity are tabulated 
in Tables IXJ - IXVIIII Statistical and systematic uncer- 
tainties are also shown. The normalization uncertainty 
(9.7% for 200 GeV, 11% for 62.4 GeV) is not included. 
For protons and antiprotons, there are two kinds of ta- 
bles, i.e. with and without the feed-down weak decay 
corrections. 
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TABLE XI: n+ and tt" cross sections [mb GeV^c^]) in p + p collisions at \/s = 200 GeV. Statistical (2nd column) 

and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (9.7%) is not 
included. 

PT [GeV/c] tt'*' 7r~ 

0.35 2.77e+01 ± 3.0e-01 ± 1.9e+00 2.63e+01 ± 3.7e-01 ± 1.8e+00 

0.45 1.45e+01 ± 1.5e-01 ± l.Oe+00 1.40e+01 ± 2.0e-01 ± 9.8e-01 

0.55 7.76e+00 ± 8.6e-02 ± 5.4e-01 7.91e+00 ± 1.2e-01 ± 5.5e-01 

0.65 4.39O+00 ± 5.3c-02 ± 3.1c-01 4.44c+00 ± 7.0o-02 ± 3.1o-01 

0.75 2.65C+00 ± 3.5o-02 ± 1.9c-01 2.69o+00 ± 4.6c-02 ± 1.9c-01 

0.85 1.59e+00 ± 2.2o-02 ± l.lo-Ol 1.60e+00 ± 2.9o-02 ± l.lo-Ol 



0.95 


l.Ole+00 ± 1.5e-02 ± 7.1e-02 


9.83e-01 


± 


1.9e-02 


± 


6.9e-02 


1.05 


6.45e-01 


± 


l.le-02 


± 


4.5e-02 


6.30e-01 


± 


1.3e-02 


± 


4.4e-02 


1.15 


4.I80-OI 


± 


7.2O-03 


± 


2.9O-02 


4.36C-01 


± 


9.5O-03 


± 


3.I0-O2 


1.25 


2.76O-01 


± 


5.O0-O3 


± 


1.9O-02 


2. 790-01 


± 


6.3O-03 


± 


2.O0-O2 


1.35 


1.88C-01 


± 


3.6C-03 


± 


1.3e-02 


I.9O0-OI 


± 


4.4O-03 


± 


1.3O-02 


1.45 


1.29e-01 


± 


2.6e-03 


± 


9.0e-03 


1.29e-01 


± 


3.1e-03 


± 


9.0e-03 


1.55 


9.07e-02 


± 


1.9e-03 


± 


6.4e-03 


9.05e-02 


± 


2.3e-03 


± 


6.3e-03 


1.65 


6.52e-02 


± 


1.4e-03 


± 


4.6e-03 


6.47e-02 


± 


1.7e-03 


± 


4.5e-03 


1.75 


4.48e-02 


± 


9.9e-04 


± 


3.1e-03 


4.69e-02 


± 


1.2e-03 


± 


3.3e-03 


1.85 


3. 450-02 


± 


8.I0-O4 


± 


2.4C-03 


3.4O0-O2 


± 


9.3O-04 


± 


2.4O-03 


1.95 


2.49e-02 


± 


6.1e-04 


± 


1.7e-03 


2.56e-02 


± 


7.4e-04 


± 


1.8e-03 


2.05 


1.83e-02 
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4.7e-04 


± 


1.3e-03 


1.81e-02 
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5.5e-04 


± 


1.3e-03 


2.15 


1.37e-02 


± 


3.8e-04 


± 


9.6e-04 


1.33e-02 


± 


4.3e-04 


± 


9.3e-04 


2.25 


1.13e-02 
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3.5e-04 


± 


7.9e-04 


1.03e-02 


± 


3.6e-04 


± 


7.2e-04 


2.35 


8.2I0-O3 
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2.80-04 


± 


5.7C-04 


7.48e-03 
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2.8e-04 
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5.2O-04 


2.45 


6.73O-03 
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2.5O-04 
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4.7O-04 


6.34e-03 
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2.7e-04 
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4.4O-04 


2.55 


5.39O-03 
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2.3O-04 
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3.80-04 
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3.02e-03 


± 


1.6e-04 


± 
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2.23O-03 
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2.95 


1.82C-03 
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1.2C-04 
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1.3e-04 


1.66e-03 
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1.2e-04 
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1.2e-04 



TABLE XIL K+ and K' cross sections {E^ [mb GeV^c^]) in p+p collisions at ^/s = 200 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (9.7%) is not 
included. 
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0.45 


1.96O-I-00 ± 


5.0C-02 ± 


1.4O-01 


1.89o-f00 ± 7.0C-02 ± 1.3o-01 


0.55 


1.35e-|-00 ± 3.0C-02 ± 


9.4O-02 


1.37O-I-00 ± 4.3C-02 
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9.60-02 


0.65 


8.71e-01 
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1.9O-02 
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6.I0-O2 


8.28C-01 
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2.3O-02 
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5.80-02 
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5.86e-01 
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1.3e-02 
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4.1e-02 


5.60e-01 
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1.6e-02 
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0.85 


3.95e-01 
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8.7e-03 
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2.8e-02 


3.87e-01 
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l.le-02 
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7.3e-03 


± 


1.8e-02 


1.05 


1.72e-01 


± 


3.9e-03 


± 


1.2e-02 


1.83e-01 


± 


5.5e-03 


± 


1.3e-02 


1.15 


1.26e-01 


± 


3.0e-03 


± 


8.9e-03 


1.16C-01 


± 


3.5O-03 


± 


8.I0-O3 


1.25 


8.52e-02 


± 


2.1e-03 


± 


6.0e-03 


8.97e-02 


± 


2.8e-03 


± 


6.3e-03 


1.35 


6.08e-02 


± 


1.5e-03 


± 


4.3e-03 


6.23e-02 


± 


2.0e-03 


± 


4.4e-03 


1.45 


4.59e-02 


± 


1.2e-03 


± 


3.2e-03 


4.27e-02 


± 


1.4e-03 


± 


3.0e-03 


1.55 


3.29e-02 


± 


9.0e-04 


± 


2.3e-03 


3.21e-02 


± 


l.le-03 


± 


2.2e-03 


1.65 


2.39C-02 


± 


6.60-04 


± 


1.7O-03 


2.23C-02 


± 


7.4O-04 


± 


1.60-03 


1.75 


1.86e-02 


± 


5.3e-04 


± 


1.3e-03 


1.81e-02 


± 


6.2e-04 


± 


1.3e-03 


1.85 


1.49e-02 


± 


4.4e-04 


± 


l.Oe-03 


1.36e-02 


± 


4.7e-04 


± 


9.5e-04 


1.95 


1.13e-02 


± 


3.5e-04 


± 


7.9e-04 


1.03e-02 


± 


3.7e-04 


± 


7.2e-04 
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TABLE XIII: p and p cross sections {E^ [mb GeV^c^]) in p + p collisions at = 200 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (9.7%) is not 
included. Feed-down weak decay corrections are not applied. 



PT [GeV/c] 






P 










p 






0.55 


1.02e+00 ± 2.0e-02 ± 


6.2e-02 


7.88e-01 


± 


1.6e-02 


± 


5.5e-02 


0.65 


7.40e-01 


± 


1.4e-02 


± 


4.5e-02 


6.04e-01 


± 


1.2e-02 


zb 


4.2e-02 


0.75 


5.58e-01 


± 


l.le-02 


± 


3.4e-02 


4.62e-01 


± 


9.1e-03 


± 


3.2e-02 


0.85 


3.77e-01 


± 


7.7e-03 


± 


2.3e-02 


3.18e-01 


± 


6.3e-03 


± 


2.2e-02 


0.95 


2.73e-01 


± 


5.9e-03 


± 


1.6e-02 


2.18e-01 


± 


4.4e-03 


± 


1.5e-02 


1.05 


1.80e-01 


± 


4.0e-03 


± 


l.le-02 


1.58e-01 


± 


3.3e-03 


± 


l.le-02 


1.15 


1.27e-01 


± 


2.9e-03 


± 


7.6e-03 


1.08e-01 


± 


2.4e-03 


± 


7.6e-03 


1.25 


9.18e-02 


± 


2.2e-03 


± 


5.5e-03 


7.54e-02 


± 


1.7e-03 


± 


5.3e-03 


1.35 


6.24C-02 


± 


1.6C-03 


± 


3.7O-03 


5.58e-02 


± 


1.3O-03 


± 


3.9C-03 


1.45 


4.80e-02 


± 


1.3e-03 


± 


2.9e-03 


3.73e-02 


± 


8.9e-04 


± 


2.6e-03 


1.55 


3.32e-02 


± 


9.1e-04 


± 


2.0e-03 


2.68e-02 


zt 


6.6e-04 


± 


1.9e-03 


1.65 


2.31e-02 


± 


6.5e-04 


± 


1.4e-03 


1.93e-02 


± 


4.9e-04 


± 


1.4e-03 


1.75 


1.70e-02 


± 


5.0e-04 


± 


l.Oe-03 


1.39e-02 


± 


3.7e-04 


± 


9.8e-04 


1.85 


1.17C-02 


± 


3.6C-04 


± 


7.0C-04 


9.69C-03 


± 


2.6e-04 


± 


6.8O-04 


1.95 


8.98C-03 


± 


2.9C-04 


± 


5.4C-04 


6.94C-03 


± 


1.9e-04 


± 


4.9O-04 


2.05 


6.68C-03 


± 


2.3C-04 


± 


4.0O-04 


5.12e-03 


± 


1.5e-04 


± 


3.6O-04 


2.15 


4.62e-03 


± 


1.6e-04 


± 


2.8e-04 


3.61e-03 


± 


l.le-04 


± 


2.5e-04 


2.25 


3.91e-03 


± 


1.5e-04 


± 


2.4e-04 


2.90e-03 


± 


9.2e-05 


± 


2.0e-04 


2.35 


2.63e-03 


± 


l.Oe-04 


± 


1.6e-04 


2.09e-03 


± 


7.0e-05 


± 


1.5e-04 


2.45 


1.79e-03 


± 


7.4e-05 


± 


l.le-04 


1.58e-03 


± 


5.5e-05 


± 


l.le-04 


2.55 


1.62e-03 


± 


7.0e-05 


± 


l.Oe-04 


l.lOe-03 


± 


4.2e-05 


± 


7.8C-05 


2.65 


1.15e-03 


± 


5.4e-05 


± 


7.2e-05 


8.85e-04 


± 


3.7e-05 


± 


6.3e-05 


2.75 


O.O^C U^i 


-1- 


4.4e-05 


-1- 


5.6e-05 


6.22e-04 


± 


2.8e-05 


± 


4.4e-05 


2.85 


6.38e-04 


± 


3.5e-05 


± 


4.1e-05 


5.07e-04 


± 


2.4e-05 


zb 


3.6e-05 


2.95 


4.97e-04 


± 


3.0e-05 


± 


3.2e-05 


3.80e-04 


± 


2.0e-05 


zb 


2.7e-05 


3.05 


4.13e-04 


± 


2.6e-05 


± 


2.7O-05 


3.13e-04 


± 


1.7e-05 


zb 


2.3O-05 


3.10 


3.80e-04 


± 


1.8e-05 


± 


2.5e-05 


2.75e-04 


± 


l.le-05 


zb 


2.0e-05 


3.30 


2.33e-04 


± 


1.3e-05 


± 


1.6e-05 


1.92e-04 


± 


9.0e-06 


zb 


1.4e-05 


3.50 


1.57e-04 


± 


l.Oe-05 


± 


l.le-05 


1.12e-04 


± 


6.5e-06 


zb 


8.6e-06 


3.70 


l.lle-04 


± 


8.9e-06 


± 


8.3e-06 


7.16e-05 


± 


5.2e-06 


zb 


5.8e-06 


3.90 


7.25C-05 


± 


7.2C-06 


± 


5.8O-06 


4.40C-05 


± 


4.0C-06 


± 


3.80-O6 


4.10 


6.23C-05 


± 


6.7C-06 


± 


5.3O-06 


3.81C-05 


± 


3.9C-06 


± 


3.60-O6 


4.30 


3.83C-05 


± 


5.5C-06 


± 


3.6C-06 


2.63C-05 


± 


3.3C-06 


± 


2.8C-06 


4.50 


3.22e-05 


± 


5.2e-06 


± 


3.3e-06 


1.82e-05 


± 


2.8e-06 


zb 


2.2e-06 
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TABLE XIV: p and p cross sections [mb GeV ^c^]) in p + p collisions at = 200 GeV. Statistical (2nd column) 

and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (9.7%) is not 
included. Feed-down weak decay corrections are applied. 



Pt [GeV/c] p p 



0.55 


5.93e-01 


± 


l.le-02 


± 


1.4e-01 


4.56e-01 


± 


9.2e-03 


± 


l.le-Ol 


0.65 


4.45e-01 


± 


8.4e-03 


± 


9.4e-02 


3.63e-01 


± 


7.0e-03 


zb 


7.8e-02 


0.75 


3.47e-01 


± 


6.9e-03 


± 


6.6e-02 


2.87e-01 


± 


5.6e-03 


± 


5.6e-02 


0.85 


2.42C-01 


± 


4.9C-03 


± 


4.2C-02 


2.04e-01 


± 


4.0e-03 


± 


3.6e-02 


0.95 


1.80e-01 


± 


3.9e-03 


± 


2.9e-02 


1.44e-01 


± 


2.9e-03 


± 


2.4e-02 


1.05 


1.22e-01 


± 


2.7e-03 


± 


1.8e-02 


1.06e-01 


± 


2.2e-03 


± 


1.6e-02 


1.15 


8.77e-02 


± 


2.0e-03 


± 


1.2e-02 


7.48e-02 


± 


1.6e-03 


± 


l.le-02 


1.25 


6.46e-02 


± 


1.6e-03 


± 


8.4e-03 


5.31e-02 


± 


1.2e-03 


± 


7.1e-03 


1.35 


4.47C-02 


± 


l.lc-03 


± 


5.5C-03 


4.00e-02 


± 


9.4O-04 


± 


5.1C-03 


1.45 


3.49e-02 


± 


9.4e-04 


± 


4.1e-03 


2.72e-02 


± 


6.5e-04 


± 


3.3e-03 


1.55 


2.45e-02 


± 


6.8e-04 


± 


2.7e-03 


1.98e-02 


± 


4.9e-04 


± 


2.3e-03 


1.65 


1.73e-02 


± 


4.9e-04 


± 


1.9e-03 


1.45e-02 


± 


3.7e-04 


± 


1.6e-03 


1.75 


1.28e-02 


± 


3.8e-04 


± 


1.3e-03 


1.06e-02 


± 


2.8e-04 


± 


1.2e-03 


1.85 


8.92C-03 


± 


2.7C-04 


± 


8.9C-04 


7.42C-03 


± 


2.O0-O4 


± 


7.9O-04 


1.95 


6.95C-03 


± 


2.2C-04 


± 


6.8C-04 


5.37C-03 


± 


1.5O-04 


± 


5.60-04 


2.05 


5.21C-03 


± 


1.8C-04 


± 


4.9C-04 


4.00C-03 


± 


1.2C-04 


± 


4.O0-O4 


2.15 


3.63e-03 


± 


1.3e-04 


± 


3.4e-04 


2.84e-03 


± 


8.7e-05 


± 


2.8e-04 


2.25 


3.10e-03 


± 


1.2e-04 


± 


2.8e-04 


2.30e-03 


± 


7.3e-05 


± 


2.2e-04 


2.35 


2.10e-03 


± 


8.2e-05 


± 


1.9e-04 


1.67e-03 


± 


5.6e-05 


± 


1.6e-04 


2.45 


1.44e-03 


± 


6.0e-05 


± 


1.3e-04 


1.27e-03 


± 


4.4e-05 


± 


1.2e-04 


2.55 


1.31e-03 


± 


5.7e-05 


± 


l.lc-04 


8.89e-04 


± 


3.4e-05 


± 


8.3C-05 


2.65 


9.31e-04 


± 


4.4e-05 


± 


8.0e-05 


7.19e-04 


± 


3.0e-05 


± 


6.6e-05 


2.75 


7.26e-04 


± 


3.6e-05 


± 


6.2e-05 


5.08e-04 


± 


2.3e-05 


± 


4.6e-05 


2.85 


5.23e-04 


1 

± 


2.9e-05 


1 

± 


4.4e-05 


A 1 ft*-. nA 


1 


2.0e-05 


1 


3.8e-05 


2.95 


4.09e-04 


± 


2.4e-05 


± 


3.4e-05 


3.13e-04 


± 


1.6e-05 


± 


2.8e-05 


3.05 


3.41e-04 


± 


2.2e-05 


± 


2.9O-05 


2.58e-04 


± 


1.4C-05 


± 


2.3O-05 


3.10 


3.14e-04 


± 


1.5e-05 


± 


2.6e-05 


2.28e-04 


± 


9.3e-06 


± 


2.0e-05 


3.30 


1.94e-04 


± 


l.le-05 


± 


1.6e-05 


1.60e-04 


± 


7.5e-06 




1.4e-05 


3.50 


1.32e-04 


± 


8.6e-06 


± 


l.le-05 


9.42e-05 


± 


5.4e-06 


± 


8.6e-06 


3.70 


9.35e-05 


± 


7.5e-06 


± 


8.2e-06 


6.03e-05 


± 


4.4e-06 


± 


5.6e-06 


3.90 


6.13e-05 


± 


6.1C-06 


± 


5.6O-06 


3.72e-05 


± 


3.4C-06 


± 


3.60-O6 


4.10 


5.28e-05 


± 


5.7C-06 


± 


5.I0-O6 


3.24O-05 


± 


3.3C-06 


± 


3.4O-06 


4.30 


3.26C-05 


± 


4.7C-06 


± 


3.3C-06 


2.23C-05 


± 


2.8C-06 


± 


2.60-O6 


4.50 


2.75e-05 


± 


4.4e-06 


± 


3.0e-06 


1.56e-05 


± 


2.4e-06 


± 


2.0e-06 
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TABLE XV: 7r+ and tt" cross sections {Ej^ [mb GeV-^c^]) in p + p collisions at \/s = 62.4 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (11%) is not 
included. 

PT [GeV/c] tt'*' 7r~ 

0.35 1.96e+01 ± 1.8e-01 ± 1.4e+00 2.08e+01 ± 1.5e-01 ± 1.3e+00 

0.45 1.07e+01 ± l.le-01 ± 7.5e-01 1.12e+01 ± 8.6e-02 ± 6.7e-01 

0.55 5.95e+00 ± 6.3e-02 ± 4.2e-01 5.94e+00 ± 4.9e-02 ± 3.6e-01 

0.65 3.38C+00 ± 3.9c-02 ± 2.4c-01 3.25c+00 ± 3.0c-02 ± 1.9c-01 

0.75 l.gic+OO ± 2.4C-02 ± 1.3c-01 1.92c+00 ± 2.0o-02 ± 1.2c-01 

0.85 1.13e+00 ± 1.6c-02 ± 7.9c-02 1.15c+00 ± 1.3c-02 ± 6.9c-02 



0.95 


6.86e-01 


± 


l.Oe-02 


± 


4.8e-02 


6.68e-01 


± 


8.4e-03 


± 


4.0e-02 


1.05 


4.30e-01 


± 


7.2e-03 


± 


3.0e-02 


4.06e-01 


± 


5.7e-03 


± 


2.4e-02 


1.15 


2.65O-01 


± 


4.9O-03 


± 


1.9C-02 


2.53C-01 


± 


4.0C-03 


± 


1.5e-02 


1.25 


1.66e-01 


± 


3.5e-03 


± 


1.2C-02 


1.60C-01 


± 


2.9C-Q3 


± 


9.6e-03 


1.35 


1.08e-01 


± 


2.6e-03 


± 


7.5e-03 


1.03O-01 


± 


2.1e-03 


± 


6.2e-03 


1.45 


7.20e-02 


± 


1.9e-03 


± 


5.0e-03 


6.74e-02 


± 


1.6e-03 


± 


4.0e-03 


1.55 


5.04e-02 


± 


1.5e-03 


± 


3.5e-03 


4.54e-02 


± 


1.2e-03 


± 


2.7e-03 


1.65 


3.48e-02 


± 


1.2e-03 


± 


2.4e-03 


3.07e-02 


± 


9.7e-04 


± 


1.8e-03 


1.75 


2.33e-02 


± 


9.5e-04 


± 


1.6e-03 


2.25e-02 


± 


8.3e-04 


± 


1.4e-03 


1.85 


1.58C-02 


± 


7.8O-04 


± 


l.lc-03 


1.55O-02 


± 


6.8O-04 


± 


9.3O-04 


1.95 


l.lle-02 


± 


6.7e-04 


± 


7.8e-04 


9.63e-03 


± 


5.2e-04 


± 


5.8e-04 


2.05 


7.13e-03 


± 


5.2e-04 


± 


5.0e-04 


7.23e-03 


± 


4.7e-04 


± 


4.3e-04 


2.15 


5.63e-03 


± 


5.0e-04 


± 


4.0e-04 


4.72e-03 


± 


3.9e-04 


± 


2.9e-04 


2.25 


4.22e-03 


± 


4.3e-04 


± 


3.0e-04 


3.32e-03 


± 


3.4e-04 


± 


2.0e-04 


2.35 


2.69e-03 


± 


3.7C-04 


± 


1.9C-04 


2.67O-03 


± 


3.3C-04 


± 


1.7C-04 


2.45 


1.96e-03 


± 


3.1C-04 


± 


1.4C-04 


1.75O-03 


± 


2.9C-04 


± 


l.lc-04 


2.55 


1.45O-03 


± 


3.3O-04 


± 


l.lc-04 


1.49O-03 


± 


2.7C-04 


± 


9.8C-05 


2.65 


9.07e-04 


± 


2.2e-04 


± 


7.0e-05 


1.07e-03 


± 


2.5e-04 


± 


7.3e-05 


2.75 


1.09e-03 


± 


3.0e-04 


± 


8.6e-05 


7.62e-04 


± 


2.5e-04 


± 


5.4e-05 


2.85 


6.48e-04 


± 


2.3e-04 


± 


5.3e-05 


5.10e-04 


± 


2.0e-04 


± 


3.7e-05 
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TABLE XVI: K+ and K' cross sections {E^^ [mb GeV^c^]) in p + p collisions at ^/s = 62.4 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (11%) is not 
included. 

PT [GeV/c] 

0.45 1.18e+00 ± 2.7e-02 ± 8.2e-02 1.06e+00 ± 1.9e-02 ± 7.4e-02 



0.55 


8.18e-01 


± 


1, 


.8e-02 


± 


5.7e-02 


7.48e-01 


± 


1.3e-02 


± 


5.2e-02 


0.65 


6.07e-01 


± 


1 


.3e-02 


± 


4.3e-02 


5.30e-01 


± 


9.6e-03 


± 


3.7e-02 


0.75 


3.72e-01 


± 


8 


.4e-03 


± 


2.6e-02 


3.43O-01 


± 


6.7O-03 


± 


2.4C-02 


0.85 


2.50e-01 


± 


6 


.lc-03 


± 


1.8e-02 


2.14C-01 


± 


4.6O-03 


± 


1.5O-02 


0.95 


1.73e-01 


± 


4 


.7c-03 


± 


1.2C-02 


1.40C-01 


± 


3.4C-03 


± 


9.8C-03 


1.05 


1.12e-01 


± 


3 


.3e-03 


± 


7.8e-03 


9.05e-02 


± 


2.5e-03 


± 


6.3e-03 


1.15 


7.94e-02 


± 


2 


.7e-03 


± 


5.6e-03 


6.17e-02 


± 


1.9e-03 


± 


4.3e-03 


1.25 


4.88C-02 


± 


1 


.9o-03 


± 


3.4O-03 


4.35C-02 


± 


1.5O-03 


± 


3.0O-03 


1.35 


3.41C-02 


± 


1 


.5o-03 


± 


2.4O-03 


2.84C-02 


± 


1.2O-03 


± 


2.0O-03 


1.45 


2.45C-02 


± 


1 


.2c-03 


± 


1.7C-03 


1.96O-02 


± 


9.2C-04 


± 


1.4C-03 


1.55 


1.63e-02 


± 


9 


.5e-04 


± 


l.le-03 


1.34e-02 


± 


7.6e-04 


± 


9.4e-04 


1.65 


1.28e-02 


± 


8, 


.Oe-04 


± 


9. le-04 


9.77e-03 


± 


6.2e-04 


± 


7.0e-04 


1.75 


9.56e-03 


± 


6 


.8e-04 


± 


7. le-04 


6.65e-03 


± 


4.8e-04 


± 


4.9e-04 


1.85 


6.34e-03 


± 


5, 


.4e-04 


± 


5.0e-04 


4.87e-03 


± 


4.2e-04 


± 


3.8e-04 


1.95 


5.28e-03 


± 


5, 


.le-04 


± 


4.4e-04 


3.45e-03 


± 


3.7e-04 


± 


2.9e-04 



TABLE XVII: p and p cross sections {E^^ [mb GeV-^c^]) in p + p collisions at ^/s = 62.4 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (11%) is not 
included. Feed-down weak decay corrections are not applied. 



PT [GeV/c] 






P 










p 






0.65 


4.63e-01 


± 


7.1e-03 


± 


4.2e-02 


3.09e-01 


± 


4.6e-03 


± 


2.2e-02 


0.75 


3.28e-01 


± 


5.4e-03 


± 


3.0C-02 


2.19e-01 


± 


3.6e-03 


± 


1.5e-02 


0.85 


2.49e-01 


± 


4.5e-03 


± 


2.2e-02 


1.59e-01 


± 


2.9e-03 


± 


l.le-02 


0.95 


1.69e-01 


± 


3.4e-03 


± 


1.5e-02 


l.lOe-01 


zt 


2.3e-03 


± 


7.7e-03 


1.05 


1.20e-01 


± 


2.7e-03 


± 


l.le-02 


7.50e-02 


± 


1.8e-03 


± 


5.3e-03 


1.15 


8.12e-02 


± 


2.1e-03 


± 


7.3e-03 


4.95e-02 


± 


1.4e-03 


± 


3.5e-03 


1.25 


5.81C-02 


± 


1.7C-03 


± 


5.2O-03 


3.32e-02 


± 


l.le-03 


± 


2.3O-03 


1.35 


3.95e-02 


± 


1.4C-03 


± 


3.6O-03 


2.37O-02 


± 


9.4O-04 


± 


1.7O-03 


1.45 


2.55e-02 


± 


9.9O-04 


± 


2.3O-03 


1.53C-02 


± 


7.1C-04 


± 


1.1O-03 


1.55 


1.84e-02 


± 


8.4e-04 


± 


1.7e-03 


1.07e-02 


± 


6.0e-04 


± 


7.5e-04 


1.65 


1.37e-02 


± 


7.2e-04 


± 


1.2e-03 


7.03e-03 


± 


4.7e-04 


± 


4.9e-04 


1.75 


9.31e-03 


± 


5.8C-04 


± 


8.4O-04 


4.49C-03 


± 


3.7C-04 


± 


3.1O-04 


1.85 


5.90e-03 


± 


4.4C-04 


± 


5.3O-04 


3.39C-03 


± 


3.4C-04 


± 


2.4O-04 


1.95 


4.02O-03 


± 


3.6C-04 


± 


3.6O-04 


2.12C-03 


± 


2.4C-04 


± 


1.5C-04 


2.05 


3.11e-03 


± 


3. le-04 


± 


2.8e-04 


1.58e-03 


± 


2.2e-04 


± 


1. le-04 


2.15 


1.99e-03 


± 


2.5e-04 


± 


1.8e-04 


1.04e-03 


± 


1.7e-04 


± 


7.3e-05 


2.25 


1.37e-03 


± 


2. le-04 


± 


1.2e-04 


6.99e-04 


± 


1.5e-04 


zt 


4.9e-05 


2.35 


8.94e-04 


± 


1.5e-04 


± 


8.0e-05 


5.90e-04 


± 


1.3e-04 


± 


4.1e-05 


2.45 


6.34C-04 


± 


1.3e-04 


± 


5.7e-05 


3.13e-04 


± 


1. le-04 


± 


2.2C-05 


2.55 


6.33e-04 


± 


1.4e-04 


± 


5.7e-05 


2.43e-04 


± 


8.3e-05 


± 


1.7e-05 


2.65 


4.56e-04 


± 


1.2e-04 


± 


4.1e-05 


1.80e-04 


± 


7.9e-05 


± 


1.3e-05 


2.75 


4.11e-04 


± 


1. le-04 


± 


3.7e-05 


1.74e-04 


± 


7.5e-05 


± 


1.2e-05 


2.85 


2.40e-04 


± 


9.5e-05 


± 


2.2e-05 


2.39e-04 


± 


9.1e-05 


± 


1.7e-05 


2.95 


1.63e-04 


± 


6.6C-05 


± 


1.5C-05 


6.57C-05 


± 


5.0C-05 


± 


4.7O-06 


3.10 


9.65e-05 


± 


3.7C-05 


± 


8.9C-06 


7.07C-05 


± 


2.8C-05 


± 


5.1e-06 


3.30 


9.05C-05 


± 


4.1C-05 


± 


8.5e-06 


4.14C-05 


± 


3.2O-05 


± 


3.1e-06 


3.50 


2.13e-05 


± 


1.9e-05 


± 


2.0e-06 


5.21e-05 


± 


3.2e-05 


± 


4.0e-06 
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TABLE XVIII: p and p cross sections (B^ [mb GeV^c^]) in p + p collisions at \/s = 62.4 GeV. Statistical (2nd column) 
and systematic (3rd column) uncertainties are shown for each particle species. The normalization uncertainty (11%) is not 
included. Feed-down weak decay corrections axe applied. 

Pt [GeV/c] p p 



0.65 


2.95e-01 


± 


4.5e-03 


± 


6.6e-02 


1.18e-01 


± 


1.8e-03 


± 


6.5e-02 


0.75 


2.38e-01 


± 


3.9e-03 


± 


3.8e-02 


1.20e-01 


± 


2.0e-03 


zb 


3.4e-02 


0.85 


1.96e-01 


± 


3.5e-03 


± 


2.5e-02 


1.05e-01 


± 


1.9e-03 


± 


1.9e-02 


0.95 


1.40C-01 


± 


2.8C-03 


± 


1.6C-02 


8.12C-02 


± 


1.7C-03 


± 


I.I0-O2 


1.05 


1.03e-01 


± 


2.3e-03 


± 


l.le-02 


5.91e-02 


± 


1.4e-03 


± 


6.6e-03 


1.15 


7.18e-02 


± 


1.9e-03 


± 


7.1e-03 


4.07e-02 


± 


l.le-03 


± 


4.0e-03 


1.25 


5.23e-02 


± 


1.5e-03 


± 


5.1e-03 


2.81e-02 


± 


9.4e-04 


± 


2.6e-03 


1.35 


3.60e-02 


± 


1.2e-03 


± 


3.4e-03 


2.04e-02 


± 


8.1e-04 


± 


1.8e-03 


1.45 


2.34C-02 


± 


9.1C-04 


± 


2.2C-03 


1.33C-02 


± 


6.2C-04 


± 


I.I0-O3 


1.55 


1.70e-02 


± 


7.7e-04 


± 


1.6e-03 


9.42e-03 


± 


5.2e-04 


± 


7.8e-04 


1.65 


1.27e-02 


± 


6.7e-04 


± 


1.2e-03 


6.19e-03 




4.1e-04 


± 


5.1e-04 


1.75 


8.67e-03 


± 


5.4e-04 


± 


8.1e-04 


3.97e-03 


± 


3.3e-04 


± 


3.2e-04 


1.85 


5.51e-03 


± 


4.1e-04 


± 


5.1e-04 


3.00e-03 


± 


3.0e-04 


± 


2.4e-04 


1.95 


3.76C-03 


± 


3.3e-04 


± 


3.5e-04 


1.88e-03 


± 


2.2e-04 


± 


1.5e-04 


2.05 


2.91O-03 


± 


2.9O-04 


± 


2.7e-04 


1.41e-03 


± 


2.O0-O4 


± 


l.lc-04 


2.15 


1.86e-03 


± 


2.4C-04 


± 


1.7O-04 


9.24O-04 


± 


1.5O-04 


± 


7.4O-05 


2.25 


1.28e-03 


± 


2.0e-04 


± 


1.2e-04 


6.21e-04 


± 


1.3e-04 


± 


5.0e-05 


2.35 


8.39e-04 


± 


1.4e-04 


± 


7.7e-05 


5.25e-04 


± 


1.2e-04 


± 


4.2e-05 


2.45 


5.95e-04 


± 


1.2e-04 


± 


5.5e-05 


2.78e-04 


± 


9.4e-05 


± 


2.2e-05 


2.55 


5.94e-04 


± 


1.3e-04 


± 


5.5e-05 


2.16e-04 


± 


7.4e-05 


± 


1.7e-05 


2.65 


4.28e-04 


± 


l.le-04 


± 


3.9e-05 


1.60e-04 


± 


7.O0-O5 


± 


1.3C-05 


2.75 


3.86e-04 


± 


l.le-04 


± 


3.6e-05 


1.55e-04 


± 


6.6e-05 


± 


1.2e-05 


2.85 


2.25e-04 


± 


8.9e-05 


± 


2.1e-05 


2.13e-04 


± 


8.1e-05 


± 


1.7e-05 


2.95 


1.54e-04 


± 


6.2e-05 


± 


1.4e-05 


5.85e-05 


± 


4.4e-05 


zb 


4.7e-06 


3.10 


9.06e-05 


± 


3.4e-05 


± 


8.5e-06 


6.30e-05 


± 


2.5e-05 


± 


5.2e-06 


3.30 


8.50e-05 


± 


3.8e-05 


± 


8.I0-O6 


3.69e-05 


± 


2.8C-05 


± 


3.I0-O6 


3.50 


2.00e-05 


± 


1.7e-05 


± 


2.0e-06 


4.64e-05 


± 


2.8e-05 


± 


4.0e-06 



